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Abstract 

Taking seriously the phenomenological indications for supersymmetry, we have made 
a detailed study of unified minimal SUSY, including many effects at the few percent 
level in a consistent fashion. We report here a general analysis of what can be studied 
without choosing a particular gauge group at the unification scale. Firstly, we find 
that the encouraging SUSY unification results of recent years do survive the challenge 
of a more complete and accurate analysis. Taking into account effects at the 5-10% 
level leads to several improvements of previous results, and allows us to sharpen our 
predictions for SUSY in the light of unification. We perform a thorough study of the 
parameter space and look for patterns to indicate SUSY predictions, so that they do 
not depend on arbitrary choices of some parameters or untested assumptions. Our 
results can be viewed as a fully constrained Minimal SUSY Standard Model (CMSSM). 
The resulting model forms a well-defined basis for comparing the physics potential of 
different facilities. Very little of the acceptable parameter space has been excluded by 
LEP or FNAL so far, but a significant fraction can be covered when these accelerators 
are upgraded. A number of initial applications to the understanding of the values of 
and mt, the SUSY spectrum, detectabihty of SUSY at LEP 11 or FNAL, BR(6 sj), 
T(Z bb), dark matter, etc., are included in a separate section that might be of more 
interest to some readers than the technical aspects of model-building. We formulate an 
approach to extracting SUSY parameters from data when superpartners are detected. 
For small tan/3 or large mt both mi/2 and mg are entirely bounded from above at 
O(lTeV) without having to use a fine-tuning constraint. 



E-mail: gkane@uinich.edu 
^E-mail: kolda@uinich.edu 
^E-mail: leszekR@uinich.edu 
''E-mail: jwells@umich.edu 



Contents 



1 Introduction 



2 Formalism] 



2.1 Basic Assumptions 



2.2 Liglit Tliresliold Corrections 



3 Bottom- Tau Yukawa Unification 



4 Electroweak Gauge Symmetry Breaking 



5 Supergravity-Based Constraints 



6 Procedure 



7 

7 



10 
13 
15 
16 



6.1 Choice of Independent Parameters! 16 



6.2 Running the RGEs 17 



7 Constraints 



20 

7.1 Limits from Experimental Searches 21 

21 



7.2 b 



7.3 Color and Charge Breaking 22 



7.4.1 Neutral LSP 



7.4.2 Neutralino Relic Abundance 



7.4 Lightest Neutrahno as the LSP and DM Candidate! 22 

22 

23 



7.4.3 Age of the Universe 24 

7.4.4 Dark Matter! 24 



8 Results 



8.1 Effect of Constraints 



j.1.3 Constraints from and on as(mz) 27 

SJA Constraints from EWSB! 27 



Ll.l 



25 

25 

Constraints from Experimental Searches 26 

26 



).1.2 Constraints from 6 — > 57 



8.1.5 Constraints from Avoiding Color Breaking 27 

8.1.6 Constraints from Neutralino LSP 28 



8.1.7 Constraints from the Age of the Universe 28 



j.1.8 Constraints from Requiring No Fine-tuning! 29 



8.2 



constrained Minimal PAramctcrS Space ( COMPASS] 29 

j.2.1 General Properties 29 



).2.2 Regions Favored by the Dark Matter Constraint 33 



8.3 Effect of the Full Effective Higgs Potential! 34 



9 Applications 



34 



9.1 Overview! 34 



9.2 Higgs Physics 35 



9.2.1 What is due to? 35 



9.2.2 Detection of the Higgs Boson 

|9.3 What Is the Origin of rtit! 

1 9.4 BR(fe ^Tt]] 

19.5 Detection of SUSY at LEP II and FNAL| 

|9.6 What If FNAL and LEP II Do Not Detect a Sparticle? . 

1 9.7 What If FNAL or LEP II Does Find a Spart"icle?| . . . . 

19.8 The T{Z bb) Partial Width| . . . . 

1 9.9 Neutrahno LSP as Dark Matter 

b.lO The SUSY Spectrum| 

|9.11 Physics Prejudice Enhancement of Part of Model Space? 



10 Summary and Comments 



1 Introduction 



There has recently been a significant amount of activity in the field of supersymmetric 
grand unification (SUSY GUT) and its possible implications for the existence of low-energy 
SUSY and for future SUSY searches. This renewed interest was primarily caused by the 
observation Q that LEP measurements of the gauge coupling constants seem to imply their 
(grand) unification in a supersymmetric theory with superpartners near the weak scale, 
reinforced by the awareness that several phenomenological outcomes were consistent with 
SUSY although they need not have been. It was shown ||l|, ^ ^, |5| that the couplings 
merge at the GUT energy scale Mx even in the simplest supersymmetric extension of 
the Standard Model, the so-called MSSM, while they badly fail to do so in the Standard 
Model (SM) alone. This remarkable fact has been interpreted by many as a strong hint 
for a SUSY GUT, especially since its main arch-rivals for expected physics beyond the SM, 
the composite and technicolor approaches, seem now even more disfavored by the precise 
measurements at LEP |^] . It has been argued that the unification of gauge couplings is also 
possible in some non-SUSY models ||^]. These models are, however, exuberantly complicated 
and lack other virtues. Besides, one should not forget that ordinary GUTs suffer from the 
hierarchy and naturalness problems which SUSY automatically cures. 

Certainly SUSY gauge coupling unification does not constitute a proof of SUSY, nor 
can it serve as a substitute for the direct discovery of a SUSY particle. On the other hand, 
it is clearly very encouraging and should not be ignored. In fact, initial simplified studies @| 
claimed that it should be possible to put stringent limits on Mx and the GUT value of the 
gauge coupling ax, as well as on the typical scale of supersymmetry breaking. Subsequently 
it was realized j^, ^, ^, ^ that additional effects, both around the electroweak and the 
GUT scale, may introduce significant modifications to the early results without, however, 
destroying supersymmetric unification. Several authors thus focussed on increasingly refined 
studies of the subtleties of gauge coupling unification [^, ^, ^, 10, 11, 13, 14]. Some |15, 
16, 1^, |l8| also considered the unification of the bottom and tau masses, which, in addition to 
the prediction of the correct value of sin^ 9^^, was regarded as a success of the early GUTs. In 
some of these studies it was argued that the mi, — m-r unification almost invariably implies 
a very heavy top quark. Many studies mentioned above typically did not address other 
important issues of the MSSM. (Some, for example, did not require correct electroweak 
gauge symmetry breaking.) Finally, some studies have adopted a more comprehensive 
approach |^, |l^, 21, 24, 25, 27|. The goal is to generate, simultaneously with 



gauge coupling unification, realistic mass spectra of the Higgs and SUSY particles. This is 
usually done in the framework of the MSSM coupled to the minimal supergravity model 
which relates many unknown quantities of the MSSM in terms of a few basic parameters at 
the GUT scale. Next, various experimental and cosmological limits can be applied to the 
resulting couplings and mass spectra. One can then examine whether all the constraints 
are consistent with each other and whether the SUSY partners have masses in the region 
of low-energy ( ^ 1 TeV) SUSY. This is the way we study the MSSM in this work. Similar 
approaches have been studied in the programs of Arnowitt and Nath in Refs. []l9|, 28, p9| 



and of Lopez, Nanopoulos, et al, in Refs [^, |2^, |31|] among others 



We want to stress that only such a comprehensive study can be regarded as relatively 
self-consistent. Considering gauge coupling unification alone neglects the contribution (at 
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2-loops) from the Yukawa couplings. It usually also assumes grossly oversimplified super- 
symmetric mass spectra. More importantly, if one wants to include also the running of the 
Yukawa couplings, one is faced with the problem of whether or not one can at the same 



time generate electroweak symmetry breaking |32], where the magnitude of the top Yukawa 
coupling is of crucial importance. Furthermore, in general one must take into account the 
running of all the Yukawa coupling of the third generation, as we do in the present study. 
In order to impose electroweak symmetry breaking properly one needs to run not only the 
Higgs mass parameters but in fact all the relevant SUSY parameters which will be specified 
below. Deriving spectra that are compatible with the gauge coupling unification and elec- 
troweak symmetry breaking can only be achieved if the whole set of relevant parameters 
is simultaneously evaluated. Further, numerical effects from 2-loops, the full 1-loop Higgs 
effective potential, etc., often significantly affect the results. 

Only after implementing this comprehensive approach are we able to reject the ranges 
of parameters that are either unphysical or experimentally excluded, while maintaining 
consistency with gauge coupling unification. Many of the detailed effects we include have 
important consequences. For example, the two typical solutions presented as a result of 
such an analysis in Ref. Q are no longer acceptable when the more complete analysis is 
done. 

Once we derive a self-consistent SUSY spectrum that follows from grand unification, 
we can compare it with the present experimental limits. Furthermore, we can study its 
implications for cosmology and derive additional bounds. Finally, we may establish what 
ranges (and properties) of the parameter space are compatible with all limits. Such ranges 
should then be focussed on in planning for experimental searches as the most "natural" — 
those expected in the Constrained Minimal Supersymmetric Standard Model (CMSSM). 

It should be emphasized that one reason it is worthwhile doing extensive work construct- 
ing SUSY models and analyzing their implications even though the full theory is not known, 
and the origin of SUSY breaking is not understood, is that the form of the Lagrangian at 
the GUT energy scale (~ 10^^ GeV) is very general and quite insensitive to our ignorance. 
The kinetic energy terms are not completely unique, but corrections are likely to be of order 
^gut/"Ipi and thus small |3^ ]. Apart from these, given the i?-parity conservation that we 
think is motivated by the stability of the proton and by cold dark matter, the superpo- 
tential we write is general, and so is the form of the soft terms [34]. Whether one arrives 
at the Lagrangian from supergravity or string theory, it has the same form [^5[ so long as 
quadratic divergences that would mix the high and low scales (i.e., terms which are not 
soft) are excluded. Thus anything one can learn about the Lagrangian by imposing physics 
constraints will be of general validity. Until superpartners are detected the information we 
have will not be sufficiently extensive to determine all parameters in it separately, of course, 
so one will have to make various simplifying assumptions. These assumptions can be tested 
in many ways as soon as superpartner masses and branching ratios are available. 

Some aspects of SUSY GUTs, most notably the GUT-scale corrections [O, 18, 3£] to 
the running of the gauge and Yukawa couplings and proton decay |19, jSj, can only be 
considered once a specific GUT model is selected. While we have no objections to most GUT 
gauge groups, for several reasons we would rather proceed by first learning what we can say 
without specifying a gauge group, and then by making a comparative study of GUT gauge 
groups. One reason is that there may be no unification group at all fs^]. In fact in many 
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string models the SM gauge group (perhaps enlarged by one or two {7(l)'s) is obtainable 
directly from strings in which case one has gauge coupling unification without an underlying 
gauge group unification. Also we are concerned about SU{5) as a unification gauge group 
because we think it would be an astonishing accident if SUSY was otherwise successful and 
also provided just the amount and kind of cold dark matter needed by cosmology, but either 
Nature did not use this dark matter or did not have it occur naturally ||3^ in the structure 
of the theory (as would have to be the case with SU{5) because i2-parity conservation has 
to be imposed by hand there). Further, while some groups |19, 3C] have shown that the 
proton decay constraint can be very important, others [37| have argued that the situation 
is not unambiguous. Thus we feel that it is important to maintain the distinction between 
the MSSM and the particular low-energy model one derives by choosing a specific GUT 
model. We will assume throughout this analysis only that the gauge couplings unify with 
sin"^ 9w{Mx) = 3/8 and that the theory remains SU{3) x SU{2) x U{1)y symmetric up to 
that unification scale. We are extending our approach to include a comparative study of 
implications of unification gauge groups (or no simple unification), and will report on this 
in a future publication. 

In Section 2A we briefly remind the reader of the basic assumptions underlying the 
MSSM. In Section we take an initial approach to the issue of gauge coupling unification 
and focus in particular on the effect of light mass thresholds. In Section |3| we digress on the 
issue of TJib — m-r unification and discuss to what extent it requires a very heavy top quark. 
The dynamical radiative electroweak symmetry breaking and the resulting constraints are 
treated in Section ^. In Section ^ we briefly list supergravity-induced relations between the 
parameters of the model, and in the next Section we use them to specify the list of indepen- 
dent parameters that we choose to perform our numerical studies. Also in this Section we 
describe the technical aspects of the procedure used in this analysis. In Section |^ we discuss 
several experimental and cosmological limits which we use in Section ^ to constrain the 
remaining parameter space. In Section |^ we also survey a number of results of our analysis 
concerning the resulting patterns of SUSY spectra. From the phenomenological point of 
view we arrive at a COnstrained Minimal PArameterS Space (COMPASS) such that every 
choice of constrained parameters is guaranteed to have gauge coupling unification, elec- 
troweak symmetry breaking, and all experimental constraints and cosmological constraints 
satisfied. COMPASS will be our guide to what predictions could really occur and are not 
excluded by any known constraint. 

COMPASS still does not uniquely determine each parameter {mi/2, rriQ, rrit, tan/?, Aq, 
sgn/xo as defined in Section ^). They can take a range of (highly correlated) values, though 
remarkably it typically implies mass spectra within the 1 TeV mass range. We want to avoid 
further assumptions about the parameters because no further theory or data is available 
to guide us, so we explore the general implications resulting from COMPASS by varying 
all relevant parameters over wide ranges of values. In future work we will explore in detail 
predictions for hadron (FNAL, LHC) and electron (LEP, LEP II, NLC) colliders, including 
to what extent SUSY is detectable at LEP II and FNAL (with upgrades) ; in Section |^ we 
give a first survey. We also study such issues as what gives the dominant contributions to 
nifi and nit, the spectrum of superpartners and predictions of SUSY for the cosmological 
abundance of the lightest supersymmetric particle (LSP), BR(6 57) and T{Z bb). In 
addition we briefly illustrate a new approach to extracting SUSY parameters from data. 
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Solving the equations giving the parameters of the Lagrangian in terms of experimental 
observables can be difficult and misleading if approximations are introduced, but with our 
CMSSM the basic parameters can be easily extracted. Section |9| can be read independently 
of the rest of the paper, and those more interested in the phenomenological implications 
rather than the technical aspects of model-building may prefer to do so. Although this 
paper is long, we think it is very important to present a single treatment that generates 
solutions of the CMSSM consistent with all theoretical and experimental constraints, and 
examines their consequences and predictions. 

2 Formalism 

2.1 Basic Assumptions 

Several features make the Minimal Supersymmetric Standard Model (MSSM) a particularly 
interesting extension of the Standard Model. The model is based on the same gauge group 
as the SM, and its particle content is the minimal one required to implement supersymmetry 
in a consistent way. It is described by the i?-parity conserving superpotential 

W = h^jQ,Huu''j + hfjQ.Hdd] + hf^L^Hae'^ + fiHaHu- (1) 

Here Q, L represent the quark and lepton SU{2) doublet superfields, u'^, (F, the corre- 
sponding SU{2) singlets, and Hu, Ha the Higgs superfields whose scalar components give 
mass to up- and down- type quarks/leptons respectively. Generational indices have been 
shown explicitly, but group indices have been dropped. In addition, one introduces all the 
allowed soft supersymmetry-breaking terms. These are given by 

+ Bii [HdHu + h.c.)+m]jJ^Hd? + m]jJ^Hu\'^ 

+ {^M^i^B^B + \M2ij'^rw + ^"^^gV-g + h.c)j . (2) 

Here the tilded fields are the scalar partners of the quark and lepton fields, while the ■(/'i are 
the spin-^ partners of the i = ?7(l)y, SU{2)l, SU{3)c gauge bosons. The ^(ij), B, and all 
other new parameters in Csoft are a priori unknown mass parameters. 

The full Lagrangian consists of the kinetic and gauge terms (which are assumed to be 
minimal), the terms derived from the superpotential (the F-terms), and Csoft- It is im- 
portant to understand that the Lagrangian we study has the most general set of i?-parity 
conserving soft-breaking terms, that is, terms that do not induce quadratic divergences 
and thereby preserve the existence of two disparate mass scales. We require i?-parity con- 
servation motivated not only by the lack of fast proton decay in nature, but also by the 
natural success of the theory in predicting the existence of dark matter. In Section ^ we 
will add some other assumptions that relate various soft-breaking terms; these assump- 
tions are somewhat motivated and can easily be removed for further study if theoretical or 
phenomenological opportunities exist. 
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The model as defined by Eqs. (||), @ is tlie simplest phenomenologically viable super- 
symmetric extension of the SM. It is also general in the sense of allowing the most general 
form [34 1 of soft terms in Eq. (^). On the other hand, because of the large number of new 
unknown parameters the model is not very predictive. A natural way of relating them is to 
think of the MSSM as coming out of some underlying GUT (or string) model. 

One possible approach is to select at the start a specific GUT which at low-energy 
would take the form of the MSSM (plus possibly a modified neutrino sector which we 
neglect here). This can be done with any GUT which can break into the SM gauge group, 
the minimal SU{5) being the simplest and most often studied choice. In this approach, 
however, one must also consider the whole GUT-scale structure with a more complicated 
Higgs sector. Guided by minimality, one often focuses on the simplest Higgs sector of 
SU{5). But that model cannot be regarded as realistic or particularly attractive due to the 
well-known problem of doublet-triplet splitting. Fixing this new "fine-tuning" problem at 
the GUT scale requires significant modifications of the model. In other words, at present 
we believe there is no commonly accepted "standard" GUT model. 

Another approach is to treat the MSSM as an effective model that could arise from a 
large class of GUT models while not making any specific choice. Instead, one can make 
various reasonable assumptions at the GUT scale consistent with general properties of that 
class of GUT models and next study "corrections" due to a specific GUT. In this approach 
one therefore initially neglects all possible corrections due to the superheavy states. This 
is the approach that we will follow here. We will be adopting more and more assumptions 
at the GUT scale, starting in the next section from just gauge coupling unification and 
eventually considering the MSSM in the framework of the minimal supergravity model. 
While we won't choose any specific GUT we will remark below about the importance of 
some of the possible corrections at the GUT energy scale. We feel it is important to 
distinguish what we can learn from this approach from the results that would be obtained 
if we chose a specific unification gauge group. 



2.2 Light Threshold Corrections 

We first address several issues that can be studied without necessarily introducing further 
simplifications of the parameter space. We begin by focussing on the running of the gauge 
couplings alone and in particular on the important role played by the mass thresholds due 
to the Higgs and supersymmetric particles. 

In running the Renormalization Group Equations (RGEs) between the weak and GUT 
scales the coefficients of the RGEs change at each particle's mass threshold due to the de- 
coupling of states at scales above their masses. Initially, a simplified case was considered 
where one assumed mass degeneracy for all the sparticles (along with the second Higgs 
doublet) at some scale usually denoted MsusY- In that case one uses the /3-functions for 
the gauge couplings of the SM between Q = mz and Q = MsusY, and those of the MSSM 
between Q = MgusY and Q = Mx- 

However, the effects of a non-degenerate SUSY spectrum on the gauge coupling /?- 
functions provides a significant correction to the naive solutions of the RGEs |, g ||, H . 
The assumption that MsusY could represent some average sparticle mass for highly non- 
degenerate spectra, such as one gets in super-unified models, is in general incorrect and can 
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lead to significant errors on the order of 10% or more in as{mz), Mx, etc. Instead one must 
take into account the various sparticle thresholds individually, changing the gauge coupling 
/^-function coefficients for each sparticle as the energy scale crosses its (running) mass, i.e., 
when it decouples from the RGEs. Accounting for each particle's contribution to the gauge 
/3-function coefficients, one can write at 1-loop [S, M: 



^MSSM ^ 4 ^ SM 2 J_ 

1 3 ^ 10 ^ 5 ^ 10 ' 



,MSSM = -11 + Iat^ + 2 ^ + 1 E {^.,^ + 0,^^ + + } (5) 

i 



where 



^ = -^a,^ + 2-loops, t = log(Q/m^), ai = |ay, 6^ = 0{Q'^ - ml). (6) 

at ZTT o 

In the summations, i = 1, . . . , A'g where Ng = 3 is the number of fermion generations, and 
= 1 is the number of SM Higgs doublets. Here also H represents the (mass degenerate) 
higgsino fields, W the partners of the IF-bosons {m^ = M2), and g the partner of the 
gluon, all taken to be mass eigenstates in this approximation. H2 is to be understood as the 
second Higgs doublet in the approximation where Hi is the SM Higgs doublet containing 
the neutral CP-even Higgs boson with mass ~ mz- H2 is heavy with each component's 
mass equal to that of the pseudoscalar Higgs. In this approximation the mixing of the two 
Higgs doublets is suppressed by inverse powers of the heavy Higgs bosons' masses and are 
therefore ignored as being of higher order and numerically negligible ||4^ . (The full 2-loop 
gauge coupling /3-functions for the SM and MSSM which we use in actual calculations can 
be found in Refs. and respectively. A discussion of 2-loop thresholds can be found 
in Section |6.2| .) 

The effect of multiple mass threshold effects on the running of the gauge couplings 
has been extensively studied recently. Notably, in a semi-analytic approach developed by 
Langacker and Polonsky [10| the effects of the thresholds on the 1-loop gauge /3-functions 
were studied. They argued that in the calculation of as{mz) from sin^ Ow, ct, and the GUT- 
unification condition, the net effect of all low-energy thresholds could be expressed in terms 
of a single scale Mg^gy (called ^SUSY in Ref. |l^). One can express this scale in terms of 
all the supersymmetric masses: 

A A 28 

M-ff (niH.Y (M2Y (M2Y ( 




' m~ mr 



9_ 

' 57 



m~ ni~ m~ 

I'eL t^iiL ^tL 



" 38 
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In the simplified case in which the spectra of squarks and sleptons are each assumed 
mass-degenerate, and taking only the contributions with leading exponents, Eq. (P) reduces 
to a similar formula given in Ref. |17|. By using a very crude parametrization in which 
rrijj — \f^\, one finds that 

^s({y ^ 1^1 ( 1 ^ 1/^1/5. (8) 

This strong dependence on fi is somewhat unexpected considering that fi does not break 
supersymmetry. 

The Mg^gY formalism is useful in providing estimates of the size of the various possible 
corrections to the running of the gauge couplings. However, it is neither accurate nor 
practical in the more comprehensive approach that we will adopt below in which the running 
of gauge couplings is simultaneously considered with the running of Yukawa couplings and 
mass parameters. Using the SM RGEs between Q = mz and Q = Mg^gy and the SUSY 
RGEs for Q > Mgygy may accurately reproduce as{mz), but it will not provide the correct 
value of ax or Mx (up to 50% errors for the latter), nor will it allow one to calculate 
correctly the ratio mb/m^(Mx). Furthermore, in this scheme 2-loop corrections to the 1- 
loop value of as{mz) derived in this method can be added only in an approximate fashion. 
These 2-loop corrections are of the same order as the 1-loop threshold corrections, and in 
fact tend to increase as{mz) by 0(10%) when included (see Table |^ in Section 3^). Thus 



we will not use the technique of an effective SUSY scale except for purposes of comparison 
in Section ^. 

In the numerical analysis that we will present later, the effect of the threshold corrections 
is automatically included separately for each contributing particle, not with a single SUSY 



threshold. We will discuss this, along with some other subtleties involved, in Section p. 2 . 

Finally, several authors have also emphasized the importance of thresholds at the GUT 
scale |l^, 11, |l^, |l^, In many models, such as minimal SU{5), these corrections 



can be sizeable. In fact, they can be comparable to the corrections coming from the non- 
degeneracy of the SUSY spectrum at the low scale (see, e.g., Ref. fl^). Consideration of 
such corrections can even be used to achieve gauge coupling unification in models where 
none seemed otherwise possible, such as non-supersymmetric SO(IO) ^. Models with 
non-minimal GUT sectors often give rise to sizable corrections that can alter low-energy 
predictions |jl[ However, consideration of these corrections can only be made after (i) 
a GUT gauge group has been chosen and (ii) the GUT Higgs sector and mass spectrum 
has been decided upon. Because we wish to study the super-unified MSSM in general, 
without reference to a particular choice of GUT gauge group or spectrum, we ignore all 
such corrections and leave them for future studies of various proposed unification schemes. 



3 Bottom- Tau Yukawa Unification 

There has been much interest recently in the issue of Yukawa coupling unification within 
the framework of SUSY. In many GUT models, including minimal SUSY-5C/(5), the down- 
type components of the lepton and quark doublets reside in the same GUT multiplets and, 
assuming a particularly simple Higgs sector, their Yukawa couplings are often equal at the 
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GUT scale. The experimentally determined ratio mh/niT- ~ 3, which decreases roughly to 
one at the GUT scale, was considered one of the early successes for GUTs. More recently 
however after the precise LEP data on gauge couplings became available, it was shown that 
the bottom-tau mass unification, while consistent with SUSY-S'[/(5), was inconsistent with 
the non-SUSY case p5[ |. 

Several groups ||l6[ [T^ , |l8|| have examined b — t mass unification more precisely in 
minimal SUSY under the assumption of gauge coupling unification. These studies have 
claimed that in order to achieve b — t mass unification one must have a top quark with 
mass very near to its IR pseudo-fixed point. That is, to a good approximation h — t mass 
unification implies [16|: 

mP°^' ~ (200 GeV) sin (3. (9) 

For top quark masses favored by LEP (130 GeV < mf"'" < 170 GeV) |§ under the assump- 
tion of a light Higgs, they find that only the small regions 1 ^ tan /? ^ 2 or tan /3 ~ 60 are 
consistent with b — t mass unification. (Here rrif"^'^ refers to the so-called pole mass of the 
top quark as opposed to the running or MS mass [47|. For a clear discussion of this point 
see Ref. [0]. We will usually speak of running masses except where we specify otherwise.) 

One is led to ask: if the top quark is found to have a mass somewhere in the LEP-favored 
region, are we absolutely forced to either very small or very large values for tan P? In order to 
answer this question one must consider how stable the stated claim is to perturbations in the 
inputs of the analysis. Such questions have been briefly considered in Refs. [0, 17, |l^, 48|. 
We find that the effects of such perturbations are often understated. 

In considering how to make the MSSM consistent with a "light" top quark (i.e. one 
with mass well below that required by Eq. (|^), we find that there are several options for 
eluding the heavy top or the extreme values of tan /3 without having to give up on 6 — r mass 
unification completely. First and foremost, it must be remembered that previous attempts 
to address this issue have suffered from a common problem: they have attempted to study 
b — T mass unification while using only a single threshold approximation for the SUSY mass 
spectrum. That is, these analyses have claimed that a non-degenerate spectrum of sparticles 
can be approximated by a single effective scale. Though this is indeed possible for a study of 
Ois{fnz) consistent with gauge coupling unification (see discussion in Section ^), no single 
threshold approximation can possibly perform the same task for b — T mass unification, 
given the dependence of the Yukawa RGEs on the gauge couplings, the presence of Yukawa 
couplings in the 2-loop RGEs, and the necessary lack of knowledge about the scale of 
unification in such an approximation. 

With this caveat in mind we now begin to explore the stability of Eq. (^) to perturbations 
in the inputs to the analysis. In this section alone we shall use the very same single threshold 
approximation about which we have just warned the reader. We do so because we are only 
interested in general numerical studies that point to possible approaches to this question, 
and because we have a consistent approach in the following sections whose results do not 
depend on the single threshold approximation for the gauge couplings. 

The one scale that we use here should not be confused with the Mg^gy introduced 
earlier, for we will choose as{mz) in this case without regard to the condition of gauge 
coupling unification. This new effective scale is in fact nothing more than the naive SUSY 
scale used in the studies of Ref. |16] and in many early SUSY studies. In displaying our 
results, we will choose this effective SUSY scale to be equal to mz', once again the exact 
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value is unimportant for our general conclusions. We also have to choose a value for the 
6-quark pole mass. In this Section, we will take the range 4.7 < m^"'*' < 5.1 GeV, which is 
the 3a bound from the recent analysis of Ref. [p9|]. 

In addition to the innate error resulting from a single threshold approximation, there 
remain other simple routes by which Eq. (|^ can be modified. We find that hy (i) allowing 
corrections to the Yukawa unification, or (ii) allowing the strong coupling constant to take 
on values near the lower end of its experimental range, one can avoid the requirement of a 
heavy top. 

The first of these routes requires one to consider corrections to the requirement that 
rrih/rrir = 1 at Mx- This is because in the interesting regions of the rrit — tan/? plane, one 
finds that in general mb/rrir < 1 at Mx- Corrections could be induced through radiative 
corrections, through effects from heavy state decoupling, through non-renormalizable oper- 
ators, or simply by the scale of 6 — r mass unification becoming displaced from the scale 
of gauge coupling unification. Without choosing any particular source, such corrections 
have been considered 18, 48]. But how large must these corrections become in order to 
significantly alter the central claim of Eq. (^)? In Fig. |l], we have shown the regions con- 
sistent with mb/m-r = 1 for bottom quark masses in the range 4.7 < mf"^'^ < 5.1 GeV and 
as{mz) = 0.120 (within the solid lines). We have also shown the region for the same range 
of bottom masses, but now with corrections to Yukawa unification of 10% (dashed lines). 
This or similar plots are most often shown in the literature as evidence for the stability of 
Eq. (I). 

Although Fig. H suggests that Eq. @ is stable to a 10% correction one might also wish to 
explore the effect of varying the strong coupling constant on the Yukawa unification. Current 
measurements of as{mz) from a variety of sources indicates that 0.110 ^ as{mz) ^ 0.130. 
Values of as{mz) in the lower half of this range in combination with a 10% uncertainty in the 
GUT relation mb/rrir = 1 significantly widen the available parameter space in the nit— tan P 
plane. (Given the analyses of Ref. IQ, such low values for as{mz) should be included in 
a careful consideration of these questions.) Such an effect is shown in Fig. |2| where we 
have taken as{mz) = 0.112. It must be emphasized that such a small value for as{mz) is 
inconsistent with the simplest SUSY-GUT unification unless we require the scale of SUSY 
masses to be 0(10 TeV). In particular, we would need a very heavy higgsino. Nonetheless, 
such a small as{mz) could come from other sources, such as heavy threshold effects or non- 
renormalizable operators. Fig. ^ clearly shows that for rrit ^ 140 GeV, all values of tan/? 
consistent with perturbative unification become allowed. Essentially, combining two 10% 
effects has eliminated the constraint among m^, mt, and tan/3. 

As we have tried to emphasize, the conclusions drawn by demanding strict b — t mass 
unification can be quite strong, yet fairly small effects due to unknowns in the analysis can 
change the results considerably. Therefore we take the following approach in the remainder 
of this paper. We will always take the r-mass as given very precisely by experiment and 
use it to determine rar{Mx)- Though the experimental uncertainty to the central value of 
4.9 GeV is larger, we will do the same for the bottom quark mass. We will not demand exact 
h — T mass unification. Because we make no specific choice of GUT group or spectrum in 
this paper, we have no mechanism otherwise for escaping the constraints imposed by 5 — r 
mass unification. Yet we also understand that corrections that will come from any eventual 
choice of GUT can, as we demonstrated above, allow a larger region of parameter space to 
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become available. When we do require exact 6 — r mass unification in our full analysis, we 
find agreement with Eq. Further, even when we do not require unification, all solutions 
generated still preserve unification to about 20%. 

Because we wish this analysis to be general and to provide insights over the entire range 
of perturbatively allowed values of tan /3 in particular, we must do without precise h — t 
mass unification. At the same time we still include everything that would otherwise follow 
from imposing this unification because solutions we find in the regions of parameter space 
consistent with Eq. (|9|) do indeed lead ioh — r mass unification. In this sense, our results are 
more general than the previously cited analyses. We think it is likely that the approximate 
unification of and rrir is telling us important physics, but we think it is premature to 
draw conclusions from it. 



4 Electroweak Gauge Symmetry Breaking 

One of the most remarkable features of the MSSM is a "built-in" mechanism for dynamical 
electroweak symmetry breaking (EWSB) [^]. The renormalization group improved super- 
symmetric Higgs potential naturally breaks SU{2) x C/(l)y — > U{l)em if the top quark 
Yukawa coupling is sufficiently large compared to the gauge couplings. As we outline below 
this will allow us to reduce the number of free parameters in the theory and express some 
GUT-scale free parameters in terms of more useful low-energy ones. 

The tree level Higgs potential can be derived from the expressions for W, Eq. (|), the 
so-called D-terms, and Csoft, Eq. (^): 

Vo = + + + h.c.) + - \HTf, (10) 



.2 

mo < 



where mf 2 = ^'h^ ~^ A*^)"^! = -^A*) the phases of the fields are chosen such that 



3 

Using the RGEs, one may define the renormalization group improved tree level Higgs 
potential, Vo{Q), at any scale Q. Vo{Q) is understood to be the tree level Vq where the fields 
and coefficients have attained a scale-dependence through their 1- or 2-loop RGEs. However, 



as was emphasized in Ref. |51], in general Vo{Q) can depend strongly on the energy scale at 
which it is evaluated. In other words, minimizing Vo{Q) at, say, Q = niz and again at some 
slightly larger Q may lead to very different values of v^, Vu, and therefore tan/? = Vu/vd- 
This behavior is due to large radiative corrections coming particularly from mass splitting 
in the t — t system. If one knew the scale Q ~ 0{mz) at which these corrections were 
small, one could safely minimize V(){Q) there. However, this scale in unknown a priori. A 
much more satisfactory solution is achieved by minimizing the full one-loop Higgs effective 
potential. The full Higgs potential can be written as 

Vmsss{Q) = Vo{Q) + AV{Q), (11) 

where (see e.g., Ref. [^ ) 

1 / s\ 

AVm-^S1.-M'ho,^--\ (12) 
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is the 1-loop contribution to Vmggs and STr f{M'^) = Ej(-l)^^ (2i + l)Trf{M'^) where M 
and j are the (field-dependent) mass and the spin of a given state, and the sum is over ah 
states in the Lagrangian. 

Electroweak symmetry-breaking can occur if the following two conditions are met: (i) 
^Higgs is bounded from below (i.e. ml + > 2|m3|), and (ii) the minimum of Vyiggs 
occurs at non-zero field configurations (i.e. mfml < m'^). It was realized early that, given 
a "large" top quark mass, EWSB could be achieved radiatively |3^. That is, despite taking 
m'jj^jm'jj^, fi'^ > at Mx, requirement (ii) above can still be satisfied. For a "large" 
mt ^ 80GeV, the running of 771-2 is dominated by ht, the top Yukawa. As the scale Q 
decreases from the GUT scale, ttt-I is driven negative while 771^ and fj? remain positive. 

Minimization of Vkiggs leads to the system of equations: 

^^"'f =ml + ml tan p + ^ml cos 2/3 + Si = (13) 
dv^ 2 

^^"'f =ml + ml cot (3 - \m\ cos 2/3 + S2 = (14) 
dvi 2 

where Si^2 = (^^^l(^^\u ^'^d all terms are implicitly Q-dependent. 



sin2/3(Q)= .2,^7 2U (15) 



- 2(ni\ I^IWJ fM2\^ ^a.,.L ijyK^ 



Solving Eqs. (pi) and (ITD one finds: 

-2mi(Q) 

^^^{Q) + ^4{Q) 

and 

2'"^(^^ - tan2 /3(Q) 
We have introduced two parameters: 

A,2{Q) ^ mi^2{Q) + ^i,2{Q) = mjj^JQ) + f^HQ) + Si,2(g). (17) 

Examining Eqs. (|l^)-(|T6|), we find that, in fact, EWSB can occur for any value of nit so 
long as nit > iTT-b- Ii^ the limit 771^ approaches ttt-^ (ignoring for now the contributions beyond 
the tree level), fj,2 approaches /if from below, but is not driven negative as in the large ttt-^ 
limit. Eq. ( [l6| ) can now only be satisfied as tan/3 approaches 1 from above. One concludes 
therefore that radiative EWSB can occur for any nit > f^b^ though small mt ( ^ my/) would 
have required tan /? ~ 1 . From a rough search of the parameter space we find that although 
the condition rTi-t ^ mw is always sufficient for EWSB (assuming appropriate values for the 
other parameters), it is also necessary in order to obtain values of tan/3 > 2. 

It would be simplest if we could always minimize VHiggs(<3) at Q = mz because we know 
from experiment the value for mz{mz) in Eq. (|l^ ) above. In minimizing Vq{Q), this would 
be dangerous. But VHiggs(Q), unlike Vb(Q), is relatively stable with respect to Q, so that 
we can choose Q = mz with confidence. 

The complete forms of Si and S2 are included in Ref. |53]. It has been emphasized |]53| , 



54 1 that the use of only the leading t — t contributions to Ay can be misleading due 
to potentially large cancellations that can occur with other terms that are not included. 
Throughout our analysis, all contributions to the complete 1-loop effective potential have 
been included. Because use of the full potential requires knowledge of the complete SUSY 
spectrum, the iterative procedure that will be outlined in Section ^ is ideally suited for 
considering this issue. 
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5 Supergravity-Based Constraints 



While the phenomenology of the MSSM is sometimes studied without referring to its GUT- 
scale origin, we want to consider in this study a highly constrained SUSY scenario with 
as many well-motivated assumptions as possible. This will of course enhance predictability 
for the ranges of parameters where SUSY may be realized. (Later we can examine what 
modifications result from relaxing assumptions.) 



As we mentioned in Section 2A, a natural and often considered approach is to couple 
the MSSM to minimal A'^ = 1 supergravity from which the following set of assumptions 
emerges: 

1. Common gaugino mass mi/2: The soft SUSY breaking gaugino mass terms are equal 
to mi/2 at Mx, 

Mi{Mx) = M2{Mx) = m~{Mx) = my^- (18) 

2. Common scalar mass tuq: The soft SUSY breaking scalar mass terms contributing to 
the squark, slepton, and Higgs masses are equal to mo at Mx, 

vr?~{Mx) = mi{Mx) = ■■■ = ml^{Mx) = mljMx) ^ ml (19) 

3. Common trilinear scalar coupling Aq: The soft trilinear SUSY breaking terms are all 
equal to Aq at Mx, 

At{Mx) = A,{Mx) = Ar{Mx) = • • • = ^o- (20) 
Through the RGEs of the MSSM, assumption ( [TsD is often expressed 

Ml = ^tan2 6l„M2 ~ O.5M2, (21) 

M2 = — m~~0.3m~, (22) 
as 9 

with Ml, M2, and m~ evaluated at the electroweak scale. One also derives rn-1/2 — I.2M2 — 
0.36m~. 

Assumptions (|l^) and (|^), in conjunction with SUSY and the gauge structure, lead to 
the following expressions for the masses of the sfermions (except for the third generation 
sfermions) at the electroweak scale (see, e.g., Ref. |5^]) 



m- =m'j + ml + bj ± "^1 cos 2/3 [r3^^'^ - Q/^^ sin^ 6'„], (23) 



where fL,R is the left (right) sfermion corresponding to an ordinary left (right) fermion, 

Tg^'^ and Qf^j^ are the third component of the weak isospin and the electric charge of 
the corresponding fermion /, and the coefficients b can be expressed as functions of the 
gauge couplings at mz and are 6 ~ 6 for squarks, ~ 0.5 for left-sleptons, and ~ 0.15 for 
right-sleptons (see, e.g., Ref. p^]). Their exact values vary somewhat with different input 
parameters. 
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While the assumptions (18), (19), and (20) derive from theoretical speculations at the 
GUT scale, we want to stress that some motivation for assuming at least the common scalar 
mass is provided by experiment. The near mass degeneracy in the — K^ system implies a 
near mass degeneracy between sl and [57|. Similarly, slepton masses have to be strongly 
degenerate from stringent bounds on /i ^ 67 [^]. It is thus sensible to generalize this 
property to all the mass terms, especially since there exists a well-motivated theoretical 
framework providing it. Alternative approaches exist |35, 58 1, though we do not consider 
them in this study. We note that for almost all topics and applications only At among the 
trilinear soft terms plays a role, so in practice we did not have to impose the condition (pO[). 
The assumptions (^), (|l9|) will be easily tested with any superpartner data. 

Many past analyses have also relied on the further assumption that Bq = Aq — mo 
[Bq = B{Mx), etc.), which follows from a restricted class of SUGRA models. As has been 
shown in Ref. |5^] , even if this relation is present at tree level in the full theory, it can be 
altered dramatically as heavy states are decoupled at Mx- We do not impose this constraint 
anywhere in the analysis. 



6 Procedure 

6.1 Choice of Independent Parameters 

After making the (SUGRA-inspired) reduction of the parameter space outlined above, we 
are left with six "fundamental" input parameters at the GUT scale: mo, mi/2, Aq, Bq, hq, 
and hfo- (In addition, we include all effects due to hi, and hr in the analysis.) However, 



not all of the parameters remain independent when we impose radiative EWSB. Eq. (16) 
allows us to eliminate fi'^{mz) as a free parameter in favor of m^, though the sign of fi is 
still free. Similarly, we can eliminate B{mz) in favor of tan (3{mz) via Eq. (p^). Finally, 
given tan/? and the RGEs we can replace hto by nit- Table || summarizes our choices. 

This "mixed" set of input parameters, mo, mi/2, Aq, /iq, tan/3, and mj, has been 
commonly used in the literature because of its technical convenience. This convenience 
becomes apparent upon inspection of the system of RGEs, in which ^ and B do not affect 
the running of any of the other parameters in the low-energy effective Lagrangian. Their 
values at the weak scale may be calculated from Eqs. (|l5|) and (jl^), and run back up to 



Inputs 


mo, mi/2, ^0, sgn^o, rut, tan/? 
rur, (mfa or 
a{mz), sin^ 6^{mz) 


Outputs 


Bo, ImoI 

(^Imx O'" ™'')' ^s{mz) 
Ml, M2, masses and mixing angles of 
gluinos, neutralinos, charginos, squarks, sleptons, 
and Higgs bosons; ^x^O' BR(6 57), etc. 



Table 1: Summary of input to and outputs from our analysis. Note that the choice between 
rrib and ^Im^ depends on whether we are testing the assumption of GUT-scale Yukawa 
unification as in Section ^ or requiring physically realistic bottom quark masses as in this 
Section and those that follow. 
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Mx in order to determine and Bq. The sign of /i is scale-independent. Note that when 
we consider tan (3 in this analysis, we will always assume tan f3{Q) = tan (3{mz) for all 
Q ~ 0{mz)] this is well-motivated by the very slow running of tan/3 and the small range 
of scales over which we consider the phenomenology of the MSSM, and so introduces only 
negligible errors. 

There is another reason for the above choice of input parameters. In some schemes, it 
is possible to determine mt as an output. We feel, however, that mt should be an input 
into any routine. Current LEP data put strong constraints on mj, and direct discovery of 
the top quark at the Tevatron may be forthcoming. Thus mt will soon serve as a relatively 
well-known input parameter. Therefore, analyses that give mt as an output will not be 
efficient in exploring the parameter space consistent with a known mt- 

There is however a certain technical difficulty associated with using the "mixed" pa- 
rametrization. Some input parameters, like the Yukawa couplings of the third generation, 
the gauge couplings, and tan /?, are known or chosen at the Z-scale. But others such as 
mo, mi and Aq are chosen at the GUT scale. Furthermore, the two scales are mixed 
in the sense that we must calculate the values of Mx and ax through the running of the 
low-energy values of the gauge couplings. This running is in turn dependent on the low- 
energy mass spectrum of the SUSY particles, which depends most heavily on the values of 
m-o and mi/2 at the GUT scale. Therefore we employ an iterative numerical procedure that 
converges on a consistent solution given all the input parameters. We discuss it below. 



6.2 Running the RGEs 

We begin our numerical procedure at the electroweak scale, which we take to be mz- This 
is an obvious choice since many experimental quantities are now available at that scale. 



At Q = mz we take as input the well- measured values of the Z mass [46|, 

= 91.187 ±0.007 GeV, (24) 
the electromagnetic coupling constant, 

aimz) = ^- , (25) 

^ ' 127.9 ±0.1 ^ ' 

and the weak mixing angle, siv? 9^{mz), in the MS scheme. (The MS value of sin^ 6'„ 
at the Z-pole is defined so that sin^ 6'w cos^ 6'w = (7ra/\/2 G_F)/m|(l — Af), where the 
radiative correction function Ar depends on both mt and m^.) The current world average 
for the weak mixing angle is sin^ 9^ = 0.2324 ± 0.0008 ± 0.0003 [|lO|, where the first error 
is due to uncertainty in the value of mt and the second error is dominated by the Higgs 
mass uncertainty. Because we take mt as a known input parameter in this analysis, the 
uncertainty due to top quark mass is replaced by a functional dependence of sin^ 0^ on 



mt li: 



sin^ = 0.2324 - 1.03 x lO"'^ 



138GeV)2 ±0.0003 (26) 



One can see the dependence of the gauge couplings on this parametrization in Table ^, where 
we have shown the values of as{mz), ax, and Mx for several values of mt, using Eq. (p6[), 
for a set of sample input parameters. It is also interesting to note that this dependence of 
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mt(GeV) 


120 


145 


170 


sin^ 6^ 


.2329 


.2322 


.2314 


asimz) 


.126 


.127 


.129 


ax 


.0414 


.0413 


.0414 


Mx/10^« GcV 


1.76 


1.94 


2.26 



Table 2: Typical effect of dependence of sin^ 9^ on for = 120, 145, 170 GeV on 
otsijnz)-, OLXi and Mx- For this table we have chosen mo = w,i/2 = 200 GeV, tan/3 = 5, 
^0 = 0, and /U > 0. 



sin^ on the top quark mass leads to a strong dependence of as(jnz) on the top mass as 
well. For constant sin^ 6'w, larger values of the top quark mass tend to decrease the value of 
asijnz) by about 3% over the allowed range of mj. However, due to the strong dependence 
of sin^ 0w on mj, the value of agimz) actually increases by about 3% over the same range. 

Besides the values of the gauge couplings at Q = m^, one also needs the Yukawa 
couplings of the third generation of quarks and leptons at mz- The mass of the r is now 
very precisely known, m-r = 1776.9 it 0.5 MeV [^0[. The mass of the 6-quark, however, 
has a larger uncertainty. Following the analyses of Refs. jlj, we take the central value 
of to be 4.9 GeV. In order to determine /i^ and h^- at Q = mz we run the 

gauge couplings a and as from their experimental values at Q = mz down to the h- and 
r-mass scales using 3- loop QCD and 2-loop QED RGEs |^]. At the mass thresholds, we 
translate |47] the experimentally measured pole masses to the MS scheme and run these 
masses back up to the Z-scale. Similarly we arrive at ht{mz) by running gauge couplings 
up to the top quark mass threshold and then running ht back down to mz- 

Now we return to a careful treatment of the threshold corrections in the running of the 
gauge couplings already mentioned in Section In the present analysis, all thresholds are 
handled as an intrinsic part of the numerical routines. Because we determine the (running) 
mass of each sparticle at the scale Q = mi{Q) anyway, we can simultaneously change the 
gauge coupling /3-function coefficients to reflect the coupling or decoupling of this particular 
state. We have already argued in Section 2.2 that the RGEs must be run at 2-loops with 
correct 1-loop thresholds, which is what we do. In Table ^ we demonstrate the importance 
of both these requirements. Notice in particular that the net effect of the 2-loop running 
is to increase as{mz) by 0(10%). Also notice that had we considered proton decay in this 
analysis, we would have found that the proton lifetime coming from dimension-6 operators 
increases when using 2-loop running instead of 1-loop by a factor of ~ 5 since Mx has 
increased by 0(50%) and the lifetime scales as M^. 

When running the gauge coupling RGEs, we follow the decoupling prescription out- 
lined in Eqs. (§)-(^)- However, there are some minor simplifications and ambiguities to 
consider ||4^. First, we decouple all higgsinos at the common scale Q = fJ.{Q), binos at 
Q = Mi(Q), winos at Q = M2{Q), and the second Higgs doublet at Q = mA{Q)- For the 
top quark, one could either choose to decouple it at its mass threshold, or simply at mz] 
numerically either procedure is essentially equivalent. One other ambiguity in the 1-loop 
RGEs arises for weak isodoublets decoupling from (32- Here, because they will always ap- 
pear in T3 = ib^ pairs in the loops, we only couple the doublet when the scale is larger than 
the heavier member of the doublet. This can be seen in Eq. (H). At 2-loops many such 
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Example 1 


Example 2 




1-loop 


2-loop 


-'"SUSY 


1-loop 


2-loop 




as{mz) 


.117 


.129 


.117 + 2- loops 


.111 


.121 


.111 + 2-loops 


ax 


.0404 


.0422 


.0456 


.0380 


.0394 


.0417 


Mx/lO'" GcV 


1.50 


2.37 


4.30 


0.79 


1.18 


1.81 



Table 3: Typical effect of dependence of as(mz), ax, and Mx on 1-loop running, 2-loop 
running, and 2-loop running with the "effective scale" of Eq. (P) for two spectra of SUSY 
particles. For Case 1, we take mo = rni/2 = 100 GeV, tan/3 = 5, mt = 145 GeV, Aq = 0, 
and n > 0. Case 2 is the same as Case 1 except mo = mi/2 = ITeV. For the two cases, 



eff 
SUSY 



13, 177 GeV respectively. Recall that Afg^gy is defined to reproduce the 1-loop 
value for as{mz)- We calculate the values of ax and Mx for Mg^gy as in Ref. [17|. 



ambiguities arise; however, the effects of individual thresholds in the 2-loop RGEs are of 
higher order and can be safely ignored. Therefore we have changed the 2-loop coefficients 
with a single threshold at Q = mi/2 above which we use MSSM 2-loop coefficients and 
below which we use those of a two-Higgs doublet SM. We have checked that dramatically 
varying the scale Q at which the 2-loop coefficients are changed from their SM values to 
the SUSY ones causes typically an 0(2%) variation in Qs(m^) (Aas(mz) ^ 0.002). Thus, 
we feel that our approximation is justified. 

It is important to reiterate that we only decouple states in the running of the gauge 
couplings. This decoupling is necessary in order to determine realistic values for as{mz)- 
However, were one to decouple states, say, from the soft mass RGEs, then one would need to 
reconsider the effective Lagrangian and matching conditions at scales below each threshold, 
where this Lagrangian, its couplings, and their RGEs would no longer be supersymmetric. 
By minimizing the 1-loop effective potential with all states included down to Q = mz, we 
effectively include the contributions from their decoupling. Therefore, in all RGEs other 
than those of the gauge couplings, we have left all states coupled down to Q = mz where 
we minimize the full 1-loop effective potential. 

Once the boundary conditions at the GUT scale have been set, we run the RGEs of 
the system in order to determine the value of a parameter at any scale Q below Mx- The 
RGEs for minimal SUSY have appeared in numerous places in the literature, including 
Refs. [43, 62]; we follow essentially the conventions of Ref. ||K 



Though various authors 
have offered semi-analytic, approximate solutions to the full set of RGEs under various 
simplifying assumptions, a full analysis of the parameter space requires that the RGEs be 
solved numerically, given the level of accuracy that we are maintaining. 

The procedure that we adopt essentially consists of repeatedly running the RGEs be- 
tween Q = mz and Q = Mx until a self-consistent solution has been isolated. 

In the first iteration for any given set of input parameters, an approximate SUSY spec- 
trum is generated. The six RGEs of the gauge and Yukawa couplings, are simultaneously 
run up first to the GUT scale using the method of Runge-Kutta. We run the gauge cou- 
plings in the SUSY-consistent DR scheme as opposed to the MS scheme, and so we impose 
the matching condition for the two schemes at Q = mz • (The net effect of the scheme 
change is less than 1% however |10, |4^.) Running up, we define Mx as that point at which 
ai{Mx) = 02 (Mx) = ax- We then set as{Mx) = ax- All scalar masses are set equal to 
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mo, all gaugino masses to 5 ^-parameters to Aq. 

The RGEs for all the twenty-six running parameters (the gauge and Yukawa couplings, 
the //-parameter, and the soft mass terms) are run back from Q = Mx down to Q = mz- For 
the gauge and Yukawa couplings, 2-loop RGEs with thresholds are used throughout, while 
1-loop RGEs are used for the SUSY soft mass parameters. Along the way, we decouple 
any particle i in the spectrum from the gauge coupling RGEs at the scale Q = mi{Q). 
As described earlier, thresholds in the 1-loop gauge coupling RGEs are used to account 
for the effects of the decoupling of the various sparticles at masses greater than mz- At 
Q = mz a value for as{mz) is found consistent with unification assumptions, and the full 
1-loop effective scalar potential is minimized in order to determine the values of iJ?{mz) and 
B{mz) that produce proper EWSB. On the next iteration when the entire set of parameters 
is again run up from Q = mz to a newly determined Mx , the parameters /i and B will also 
run, providing their corresponding values at the GUT scale. 

This entire procedure is repeated several times, terminating only after changes in the 
solutions to the RGEs are small compared to the values themselves or to the experimental 
errors, whichever are relevant. Each iteration provides a more precise spectrum of sparticles, 
which in turn provides more precise running of the gauge and Yukawa couplings. We find 
that the whole procedure is extremely stable, usually converging to a solution in just a few 
iterations. 

In Figure ^, we give an example of the running of various sparticle masses from the 
GUT scale down to the electroweak scale. Notice that the mass of the Higgs that couples 
to the top quark is driven imaginary (i.e., its mass-squared is driven negative) at scales of 
0(1 TeV), signally the onset of EWSB. This is shown in the plot as the mass itself going 
"negative" for convenience of presentation. 

When the program has isolated a solution, we have as our output all sparticle masses 
and mixings valid to 1-loop, Higgs masses which include all third generation contributions 
to the 1-loop radiative corrections [S3|, as{mz)-, ax, and Mx valid to 2-loops, and the 
GUT-scale parameters Bq and /xq. 



7 Constraints 

In applying the numerical procedure described in the previous section we have required the 
gauge coupling to unify, and from the input values of q, Eq. (25), and sin^ 6'„, Eq. (|26|), ob- 



tained a range of as{mz) as a function of independent parameters. We have also demanded 
proper EWSB yielding the experimentally measured value of mz- We have parametrized 
the many mass parameters of the MSSM in the usual way, assuming common gaugino 
and scalar masses and the A parameters, Eqs. (|T8|), (p!s|), and (pO|), as implied by minimal 
SUGRA. Before we present our results in the next section, we now list and briefly elaborate 
on several other constraints that we will impose on the output of our numerical analysis. 
As we explained in Section ^, we do not impose the condition = m-j- at the GUT scale 
because the resulting bottom quark mass is likely to be very sensitive to the threshold cor- 
rections at Mx, which we cannot include without selecting a specific GUT model. Without 
such corrections we obtain the values of m^"^'^ about 20% above the current experimental 
range, except for very large mt- 
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7.1 Limits from Experimental Searches 

LEP experiments have placed lower limits on the chargino mass of about 47GeV, and on 
the charged slepton, sneutrino, and squark masses of about 43GeV [|^]. The hghtest stop 
mass bound is dependent on the left-right stop mixing, which can reduce its coupling to 
the Z-boson. DELPHI [^5| has excluded below about 45 GeV, except for a rather tiny 
range of the mixing angle which allows m^^ > 37 GeV. 

Placing an experimental lower limit on the masses of the Higgs bosons is in principle more 
complicated since either h or A, or both, can be light, and because of potentially sizable 
radiative corrections to their masses due to the heavy top quark. Assuming reasonable 
ranges of value for mt, rtif;, and tan /3 > 1 the bounds ruh > 44 GeV and ruA > 21 GeV have 



been derived by ALEPH |64]. Other LEP experiments obtained similar limits. However, 
once we impose the unification and EWSB conditions, we find that h couplings are very 
SM-like (sin2(/3 - a) 1) so that in practice the LEP limit of about 62 GeV |6|] for the 
SM Higgs applies to h as well. For related reasons, A is always heavier than mz for us, so 
that the LEP limits on A place no serious bound. 

Lower mass bounds on the squarks and gluino have been reported at 126 GeV and 



141 GeV [S7|, respectively, assuming no cascade decays. By including cascade decays one 



can reduce those bounds by some 20 GeV or more The squark masses could become 
even as light as allowed by LEP if m~ becomes large. All squark and gluino bounds are 
very model dependent. 

LEP experiments alone cannot place a lower bound on the mass of the lightest neutralino 
because its coupling to the Z can be strongly suppressed and it is not directly detectable. 
It is only by combining LEP direct chargino and neutralino searches with indirect (Z-line 
shape) searches and with the lower bound on m~ from the Tevatron (via Eq. (|22|)) that a 
bound > 18 GeV [^] can be derived for any tan/3 > 1. 

7.2 sj 

Recently, CLEO has reported an upper bound on BR(6 57) < 5.4 x 10"^ |70|. A 



central value (3.5 x 10 ^) and a lower limit (1.5 x 10 ^) is obtained from the detection of 
B K*-f |70| and assuming that the ratio of BR{B K*j) to BR{b sj) is 15% |7T|. 
In addition to the SM contribution, SUSY allows for one-loop diagrams with the exchange 
of the charged Higgs and the charginos and neutralinos |72, ItI, 74, 75, We calculate 



BR(6 — > 57) with the formulas of Ref. |73]. These use QCD corrections that are less accurate 



than has been done for the SM case recently |77|. We are in the process of combining our 
results with those of Ref. |77] to obtain improved QCD corrected CMSSM predictions. 

In Section ^ we will apply the upper bound BR(6 — > 57) < 5.4 x 10~^. In Section 
we will present the predictions of CMSSM for BR(6 — > 57) and show that the range favored 
by CMSSM naturally falls into the range resulting from the CLEO analysis. We will also 
show that the claims |Q of a stringent bound on the charged Higgs mass are too strong in 
the CMSSM. 
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7.3 Color and Charge Breaking 



In the MSSM the Higgs potential automatically conserves color and charge but the same is 
not necessarily true with the full scalar potential. If one wishes to determine the form of the 
global minima, one must numerically search for all local minima of the full scalar potential, 
including the charged and colored states, and determine the broken symmetries associated 
with each. This is outside the realm of the study we are reporting on here and so we only 
demand that the (mass)^ of any charged or colored mass eigenstate remains positive. In 
fact, this will be an important constraint in some regions of the parameter space (especially 
for large Aq) where the lighter stop (mass)^ can become negative due to a large — tn 
mass splitting. 

It is sometimes stated in the literature that a necessary condition for avoiding color- 
breaking is to demand that I^oI/^t-o < 3. However, as pointed out by Ref. |79|, this 
condition is really neither sufficient nor necessary. Therefore, we consider values for |^o|/"io 
slightly larger than three. Knowing that minimization of the full scalar potential may lead 
to color- or charge-breaking minima for such large ratios, the constraints coming from our 
analysis can only be strengthened by a full treatment of this color /charge-breaking. 

Nonetheless, as will be discussed in the next Section, never goes negative for smaller 

values of (|^o|/"io ^ 1) simply because the — mass splitting is dominated by ^o- 



7.4 Lightest Neutrahno as the LSP and DM Candidate 

In the absence of -R-parity breaking the LSP remains absolutely stable. Depending on its 
nature it may have to face potentially tight cosmological constraints which we will discuss 
below. It will also have important experimental consequences for possible SUSY signatures. 
In the MSSM, any of the superpartners could in principle be the LSP because their masses 
are virtually unrelated. In the CMSSM the picture is very different: the masses of the 
superpartners are highly correlated. These relations are determined by the assumptions of 
Eqs. (p!8[)-(pO|) and lead to a hierarchy among the sparticle masses. As a result there are 
very few possible candidates for the LSP. Typically it is the lightest of the four neutralinos 
that comes out to be the LSP, and it has been usually favored in most phenomenological and 
cosmological studies. However, for some combinations of parameters some other sparticle, 
like the stop, the stau, or the sneutrino can be the LSP. Each of the resulting types of the 
LSP must meet cosmological constraints. 

As we have already discussed in Section |7.3| , due to a large mass splitting in the ti — tR 
sector, the lighter stop mass eigenstate (ti) may in certain cases become very light. In 
fact, one may even encounter < 0. On the other hand, the lighter stau sometimes 
becomes the LSP. As concerns the sneutrino, after we apply experimental limits and reject 
unphysical cases, we never find it to be the LSP. 



7.4.1 Neutral LSP 

It would be difficult to imagine that an electrically charged or colored massive stable particle, 
like the stau or the stop, could exist in any meaningful amount in the Universe If 
it did, it would interact with photons and become detectable. It would also interact with 
ordinary matter and dissipate its energy thus falling towards the cores of galaxies. It would 
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form stable isotopes of chemical elements. For these, and other, reasons, only electrically 
neutral and colorless particles are believed to be able to exist in the Universe in the form 
of dark matter |8C, We will therefore reject those regions of the parameter space where 



either the stop or stau are the LSP. In the rest of the study we will only deal with the 
neutralino as the LSP. 

7.4.2 Neutralino Relic Abundance 

Any stable (or meta-stable) species predicted by theory would contribute to the total mass- 
energy of the Universe. A relic abundance is usually expressed as the ratio of the particle's 
relic density to the critical density pcrit = SHq/SttG = 1.9 x 10^'^^ {hQ)g/cm'^ 

Px 



n^^f^ (27) 

Pcrit 

where pcrit corresponds to the flat Universe and /iq is the present value of the Hubble 
parameter Hq in units 100 km/s/Mpc (/iq = iqq kmfs/Mpc ) • Current estimates only require 



o.4</io< 1 ig|. 



A super symmetric LSP, being stable, cannot decay on its own but can pair-annihilate 
into ordinary matter. Its relic abundance ^x^'o inversely proportional to the LSP annihila- 
tion cross section and thus depends on the masses and couplings of the final and exchanged 
particles. In calculating the neutralino relic abundance we include all the relevant LSP 
pair-annihilation channels into ordinary matter that are kinematically allowed. Lighter 
annihilate only (except for rare radiative processes) into pairs of ordinary fermions via the 
exchange of the Z and the Higgs bosons, and the respective sfermions. (We do not include 
final state gluons since the relevant cross section has been shown to be relatively insignifi- 
cant in calculating the relic abundance in the early Universe fS^.) As grows new final 
states open up: pairs of Higgs bosons, gauge and Higgs bosons, ZZ and WW, and tt, all of 
which we include in our analysis. The actual procedure of calculating the relic abundance is 



quite involved and has been adequately described elsewhere (see, e.g., Refs. [|0|, 83, Q). We 
use the technique developed in Ref. |l83[| which allows for a reliable (except near poles and 
thresholds) computation of the thermally averaged annihilation cross section in the non- 
relativistic limit and integration of the Boltzmann equation. This technique is applicable 
to calculating the relic abundance in most of the parameter space. 

As was first pointed out in Ref. |85[, and rediscovered and elaborated by Griest and 



Seckel [36|, special care must be applied to calculating the relic abundance near the poles 



of exchanged particles and when new mass thresholds become kinematically accessible. In 



particular, proper treatment of narrow poles has been provided in Refs. |23, 28, ^] and it 
was shown that standard techniques may lead to errors reaching even two or three orders of 
magnitude in the vicinity of a pole. This is especially true for the lightest Higgs because the 
width of h is extremely narrow, and also near the Z-boson pole where the effective coupling 
is somewhat stronger. We find that the regions of the parameter space where our (standard) 
calculation fails are relatively small albeit non-negligible. In presenting our results in the 
next Section we will therefore point out those regions where the presented results for the 



neutralino relic abundance are not trustworthy. It has been argued in Refs. [p3| , 28 1 that the 
regions where the h and Z-poles dominate are favored by current limits on the proton decay 
in the SUSY SU (5) model. Since we do not select SU (5) as a GUT symmetry, nor view it 
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as particularly attractive, at this point we choose not to pay special attention to calculating 
the relic abundance near the poles. We will comment on these effects in discussing results. 



7.4.3 Age of the Universe 

In the Standard Cosmological Model the age of the Universe depends on the total relic 
abundance r^TOT- Conversely, estimates of the Universe's age place a constraint on < 
r^TOT- A conservative assumption that the Universe is at least 10 billion years old (and 



ho > 0.4) leads to |0| 



^xhl ^ 1- (28) 

If the age of the Universe is at least 15 billion years, as many currently believe, then the 
bound (|2^ ) becomes much stronger: ^x^o ~ ^.25 [|^. This is because an older Universe 
corresponds to a smaller expansion rate /iq- No stable particle can contribute to f^TOT 
more than is allowed by at least the bound ( p8D without distorting the Universe's evolution. 
This constraint is independent of the nature (or even existence) of dark matter (DM) in the 
Universe. The bound of Eq. (28) must be satisfied for any choice of free parameters and, as 



we will see in the next Section, it provides a very strong constraint on the parameter space. 
7.4.4 Dark Matter 

The visible matter in the Universe accounts for about 1% of the critical density. There 
is at present abundant evidence for the existence of significant amounts of dark matter in 
galactic halos ~ 0.1) and in clusters of galaxies (0, > 0.2) [^]. Big Bang nucleosynthesis 
(BBN) constrains the allowed range of baryonic matter in the Universe to the range 0.02 < 
i^B < 0.11 (and more recently « 0.05; see the second paper of Ref. [Q). The 
value r^TOT = 1 is strongly preferred by theory since it is predicted by the models of cosmic 
inflation and is the only stable value for Friedmann- Robertson- Walker models. Values of 
r^TOT larger than those "directly" observed are also strongly supported by most models of 
large structure formation. This, along with estimates given above, implies that: (i) most 
baryonic matter in the Universe is invisible to us, and (ii) already in halos of galaxies one 
might need a substantial amount of non-baryonic DM. If J^tot = 1 then most (about 95%) 
of the matter in the Universe is non-baryonic and dark. Current estimates of ho give, for 
^TOT = 1, 0.5 ^ ho < 0.7 (the upper bound coming from assuming the age of the Universe 
above 10 billion years), in which case one expects 0.25 < r^TOT^o ~ ^.5. While it is not 
unlikely that the galactic halos consist to a large degree of various extended MACHO-type 
objects[| (like Jupiters, brown dwarfs, etc.), it would be very hard to believe that such 
objects could fill out the whole Universe without condensing into galaxies. This, along with 
the bound on baryonic matter provided by BBN has led to a widely accepted hypothesis that 
the bulk of DM in the Universe consists of some kind of weakly interacting massive particles 
(WIMPs). The relic abundance of the lightest neutralino x (most naturally of bino-type [pT| ) 
often comes out to be in the desired range thus making it one of the best candidates for 



DM 1 92]. Being nonrelativistic, it falls into the category of cold DM (CDM) which has been 



""^ Recently, a few candidate events for MACHO's with mass ~ O.IMq have been reported by microlensing 
experiments jsoj thus implying that some sort of small stars comprise a significant component of the halo 
of our Galaxy. We note that, with the present efficiency, this discovery does not, and will not for the next 
several years, be able to eliminate other kinds of candidates for the dark matter WjA. 
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favored by models of large structure formation, in contrast to hot DM (HDM), like light 
neutrinos. In a purely CDM scenario one assumes that the LSP dominates the mass of the 
Universe, leading roughly to 

0.25 < Q^hl < 0.5. (CDM) (29) 

Motivated by the theoretical expectation that SUSY GUT theories will also have massive 
neutrinos, and phenomenologically by the result that (in the aftermath of COBE) a mixed 
CDM+HDM picture (MDM) seems to fit the astrophysical data better than the pure 
CDM model, we also consider a smaller value of ^ly.. In the mixed scenario one assumes 
about 30% of HDM (like light neutrinos with niu — 6 eV) and about 65% of CDM (bino-like 
neutralino), with baryons contributing the remaining 5%. In this case the favored range for 
O^/ig is approximately given by 

0.16 < n^hl < 0.33. (MDM) (30) 

(Strictly speaking, in the MSSM the neutrinos are massless and as such could not consti- 
tute interesting HDM. But it is straightforward to extend the model to include right handed 
neutrinos (and their sneutrino partners) and give them mass terms. We don't expect this 
extension to sizably modify the running of all the other parameters of the MSSM. It is 
with this implicit assumption that we will apply the range given by (|30| ) in analyzing the 
resulting implications for SUSY searches.) 

Both scenarios assume a significant amount of LSP DM. The sneutrino, an early candi- 
date [^] for DM, is now strongly disfavored. After the LEP experiments have placed a limit 
on its mass m~ > 43 GeV, its relic abundance can now only be negligibly small (i7~ 10~'^). 
We thus find it remarkable that we never find the sneutrino to be the LSP. Had it been the 
LSP instead of the neutralino in most of the parameter space then the CMSSM would not 
have provided a viable candidate for the DM problem! 

8 Results 

We now proceed to discuss the numerical results obtained by using the procedure for gener- 
ating low-energy output described in Section ^. We will first analyze the impact of several 
experimental, theoretical, and cosmological constraints on the parameter space. Next, we 
will focus on the region of the model's parameter space consistent with all the adopted con- 
straints and discuss the resulting consequences for the value of as{mz), the mass spectra 
of the Higgs and supersymmetric particles, and other predictions. 

8.1 Effect of Constraints 

We have generated a large set of solutions for a broad range of input parameters. We 
explore wide ranges of both mi/2 and mo, each between 50 GeV and approximately 3TeV 
in 22 logarithmic steps, for discrete values of mt = 120, 145, 170 GeV, tan/3 = 1.1, 1.5, 3, 
5, 10, 15, 20, 30, 40, 50, and Aq/uiq between -3.5 and 3.5 in increments of 0.5. We also 
consider both signs of fiQ. 
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The choice of a logarithmic scale for is technically motivated. We are 

interested most particularly in lower values of the soft masses where the fine-tuning reintro- 
duced by SUSY breaking is smallest and where we can expect currently planned facilities to 
best probe the parameter space. Likewise, the difference between tan/3 = 1.5 and tan/3 = 3 
is more significant than that between tan /3 = 40 and tan /3 = 50. Lastly, the scaling of 
Aq with mo is motivated by SUGRA, with bounds motivated by the fear of color-breaking 
global minima for large Aq. 

Mass scales above 1 TeV may seem unnatural but we also wish to explore the asymptotic 
behavior of our results. For the top mass, the three representative values, mt = 120, 145, 
170 GeV, help us sample the whole region of top mass preferred by the analysis of the LEP 
data. We pay particular attention to the middle value, mt = 145 GeV, as being favored by 
LEP (when one includes a light Higgs as required by SUSY) and perhaps by the cross section 
for candidates from the Tevatron. For tan (3 we sample a spectrum of values over the entire 
region consistent with the Yukawa couplings of the third generation remaining perturbative 
all the way up to the scale of unification. Values of tan/3 < 1.1 become difficult to study 
due to a dangerous cancellation in Eq. ( |T^ ) and because we are close to the perturbative 
limit of the top Yukawa; values of tan /? ^ 50 become difficult because the bottom and r 
Yukawa couplings are likewise close to their perturbative limits. 

Overall, we explore well over 100,000 combinations, each representing a unique point 
in the space of mt, tan/3, mi/2, mQ, Aq, and sgn/XQ. We present some representative 
solutions in Figs. in the plane {mi/2,mo). In this Section we focus mostly on the case 
mt = 145 GeV and several representative choices of tan/3 and ^o- We will also display the 
dependence on mt below. 



8.1.1 Constraints from Experimental Searches 

As we can see from Figs. at present the regions of the plane {mi/2,n^o) excluded by 
direct and indirect searches for SUSY at LEP and FNAL (see Section ^[^) are limited at 
best. The strongest direct constraints on mi/2 come from m ± > 47 GeV and/or the CDF 
gluino mass bound. Assuming m~ > 141 GeV, and neglecting cascade decays, corresponds 
roughly to mi/2 ^ 50 GeV. On the other hand, in general there is no lower bound on rn-o 
except for small mi/2 from the experimental lower bounds on the slepton and squark masses. 
As an example, we present in Fig. |l^ the region of the plane {mi/2,mo) ruled out by the 
LEP bound > 43 GeV for two extreme values of tan (3 and by m~ > 141 GeV. For 
mi/2 ^ rriQ the exact value of mo becomes unimportant, as mi/2 will come to dominate 
the values of all masses and will dictate how EWSB occurs. Though for the major portions 
of this study we have taken mg > 50 GeV, we have explored the regions of much lower mo 
and found nothing to change our conclusions as reported in Sections |8| and ^. 

In addition, we find some regions where the lighter stop mass becomes smaller than 
the current experimental bound of about 37 GeV and quickly becomes tachyonic as will be 
discussed below. 



8.1.2 Constraints from h 



Following Section 7.2, we apply the upper bound BR(& 57) < 5.4 x 10 ^. Interestingly, 
this bound is often quite important and particularly probes regions of small to moderate 
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mii2 and mo (Figs. Q-^). As ^'^'^ "'-o grow, BR(6 57) tends to decrease and produce 
the range of values consistent with CLEO for a wide range of parameters as will be shown 
in Section |9.4| . 



8.1.3 Constraints from and on ^^(m^) 



As we can see from Figs. |^|^, the values of as{mz) resulting from our analysis generally fall 
into the experimentally allowed range. LEP event shape measurements alone give as{mz) = 
0.123 lb 0.006 [H6[ while other LEP analyses and low-energy experiments typically yield 



somewhat lower ranges leading to the world-average as{mz) = 0.120 it 0.006 it 0.002 |46|. 
(We note, however, that much smaller values of as{mz) = 0.107 it 0.003 have been derived 
in Ref. |^0|.) We find that as{mz) generally decreases with growing mi/2 "^Oj and 
increases with mt (see Table |2| and Section |6.2D . Since small 1121/2 ^0 are excluded by 
some experimental constraints (Section ^.1.1| ), we find asimz) ^ 0.133, including the range 
of very small rriQ. This is a significant constraint on the entire picture and an important 
prediction. No interesting upper bound on the plane (mi/2,"Zo) can be derived from a 
lower bound on as{mz) because as{mz) decreases very slowly and reaches 0.110 for mi/2 
and/or mo in the range of tens of TeV. Keeping SUSY masses below about 1 TeV provides 
a lower bound Qs(m^) > 0.119, while requiring no fine-tuning (/ < 50, see Section |8.1. 



gives as{mz) ^ 0.118. Clearly, as the graphs show, larger values of as{mz) are favored by 
low-energy SUSY. Finally, if as{mz) comes out much smaller (as claimed in Ref. [50|), one 
may have to necessarily include GUT-scale corrections to the running of the gauge couplings 
in order to possibly remain in the theoretically favored region of low-energy SUSY below 
a few TeV. This may actually come out to be a useful tool in an attempt to discriminate 
between different GUT scenarios. 



8.1.4 Constraints from EWSB 

Proper EWSB is not automatic and requiring it places additional strong constraints on the 
allowed combinations of parameters. As can be seen in Figs. this constraint excludes 
significant regions in the upper left-corner (mo ^ "^1/2) of the plane imif2,rno), unless 
tan P is close to one or ^o is larger and negative. For > there are additional regions in 
the lower right-hand corner (mi/2 ^ 1^0) of the plane (mi/25"^o) which are also excluded 
for larger values of tan/3. This is because the full 1-loop effective potential has become 
unbounded from below in those regions. 



8.1.5 Constraints from Avoiding Color Breaking 

As we said above, sometimes m?^ becomes negative. As one can see from the presented 

figures, this usually happens roughly for mo ^ rni/2 for rather large values of |^o| (see 
symbol "L" in these areas in, e.g., Figs. ^c-d). The regions where m?^ < always grow with 

increasing tan/?. More specifically, for /io < 0, color breaking occurs when Ao/niQ < — 2 
for the whole range of tan/3, and also to some extent for Aq/itio > 3 and large tan/3. For 
the smaller values of Aq is always positive, as expected. For /io > the situation is 
generally similar for a reversed sign of ^o- 
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8.1.6 Constraints from Neutralino LSP 



As we have argued in Section 7.4.1 , only the lightest neutrahno LSP remains a viable 
candidate for DM. On the other hand, for mx/2 ^ "^o we invariably find that the hghter 
stau is the LSP, and not the neutralino, as one can see in Figs. This is expected since 
the mass of the neutralino is given roughly by ~ Mi ~ 0.4m^/2- the other 
hand, the mass of the hghter stau tr (see Eq. (^)) grows somewhat more slowly with 
™-i/2) ~ 0.38m]^/2- III region of large mi/2 ( ^ 400 GeV) and small mo, tr (and in 
fact also cr and Jir) become lighter than the lightest neutralino. For a fixed mi/2, as mo 
grows, so does m~ and x becomes the LSP again. Insisting on the neutralino LSP provides 
a very important constraint on the plane {mi/ 2, "mo), excluding the region mi/2 ^ i^o- We 
note, however, that the regions where x is LSP correspond to large m~ > 1 TeV. 

Also, we never find the sneutrino to be the LSP: regions of small mi/2 where this could 
take place have been excluded by LEP. We thus find that, in the most interesting region 
of low-energy SUSY it is the neutralino which is most often the LSP. It is also mostly 
gaugino-type (bino-type) - this will be discussed in more detail in Section ^.2| . 

8.1.7 Constraints from the Age of the Universe 

For gaugino-like x's the relic abundance ^x^o depends most strongly on the mass of the 
lightest exchanged sfermion in xx ~^ ff'-i roughly ^x^o m^f/rn^ [^4|. All sfermion masses 
grow with increasing mo, and in the case of sleptons much more slowly with mi/2, so one 
expects that the bound ^x^o ~ ^ which results from requiring that the age of the Universe 



be at least 10 billion years (see Section 7.4.3), will be stronger for mo than for mi/2- This is 
indeed often the case in the remaining regions of the parameter space. The constraint ( ^8|) 
excludes large values of mo roughly above 1 TeV and often even above a few hundred GeV. 

For small tan/3 (tan/3 w 1), the bound ( pSj ) is typically much stronger and excludes 
mo ^ 300 GeV and mi/2 ~ ITeV. As tan/3 grows slightly to at least moderate values (2 
and above), the bound becomes less constraining primarily for Aq around zero or positive 
allowing for somewhat larger values of mo and also opening the region mo ~ ^11-1/2 above 
1 TeV. This is because the s-channel Z-exchange in the process xX ~^ ff and the x pair 
annihilation into pairs of of light Higgs bosons h become unsuppressed and can reduce the 
LSP relic abundance. The Z pole effect is clearly visible in the region mo ^ 1TI1/2 — 120 GeV 
(see, e.g., Figs. ^d, or ^). But it is also in the region near this pole (and likewise near 
the h pole) that the exact calculation of the relic abundance becomes difficult. We have 



highlighted these regions in Fig. 11 



The process xX ~^ is rarely dominant but it can reduce the relic abundance con- 
siderably, especially in the most interesting region of mx/2 and mo in the range of a few 
hundred GeV for larger values of tan /3. This is clearly visible in Fig. ^ (see also Fig. |l5|) 
where the region to the right of an "island" of ^x^o ^ ^ (large mo and mi/2 ~ 270 GeV) is 
again allowed because the final state hh becomes kinematically allowed. This effect is not 
present for small tan /3 1 (compare Fig. |^ because the coupling hxx vanishes there. 

Overall, the bound ^x^o ~ ^ typically provides a very stringent constraint on the regions 
of the parameter space not already excluded by other criteria. It excludes mo roughly above 
1 TeV, except for large m^/2 where some SUSY sparticle masses (e.g., m~) become very much 
larger than 1 TeV and are therefore disfavored by the fine-tuning criterion. 
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8.1.8 Constraints from Requiring No Fine-tuning 

Finally, it is clear that if SUSY is to replace the SM as an effective theory at the electroweak 
scale, its mass parameters should not be much larger than mz- Stated differently, since 
the combination of and in Eq. (H) has to give m^, one would have to tune those 



parameters to a high precision, unless they were broadly within a ITeV mass range [95|. 
This fine-tuning in the potential minimization is a remnant of the fine-tuning exhibited 
by the full theory. In the full theory, one would parametrize fine-tuning most naturally 
by / = A|.^^y/m|. Instead, because radiative EWSB connects the SUSY scale to the 
electroweak scale, we choose to parametrize it by 

f^\m\\/ml (31) 

which is particularly stable in terms of the running of the RGEs and the minimization 
of the 1-loop effective Higgs potential. (At the tree-level our definition is similar but not 
identical to the definition of Ross and Roberts |P.) The concept of fine-tuning is somewhat 
subjective and various authors have used different definitions and criteria. 

Fig. |l^ shows the typical scaling of the fine-tuning constant with the scale of SUSY for 
a sample choice of input parameters. In order to exclude regions where large fine-tuning 
must be invoked, we will later place an upper bound of / < 50. As we can see from 



Fig. 13, this criterion typically selects the heaviest sparticle masses below roughly 1 TeV. 



It is worth stressing however that, for large tan /3, both m~ and m~ can be significantly 
larger without any excessive fine-tuning. Thus simple cuts m~, m~ < 1 TeV often made in 



the literature [ p!9| , |2l| , p2| may in general be too strong. 

One might hope that physics constraints would eliminate the need for adding a separate 
fine-tuning constraint. That indeed is the case for large ranges of parameters, which is very 
encouraging. For example, for large rrit = 170 GeV we find that the constraint ^x^o ^ 
cannot be satisfied if tuq or mi/2 are larger than several hundred GeV. This is also true for 
smaller mt if tan/3 is close to one. In general much larger mo and mi/2 become allowed as 
tan /3 grows, but this does demonstrate the kind of argument that might lead to physical 



constraints on the parameter space in place of fine-tuning |20|. 

We will not apply the constraint / < 50 in the rest of this section because we also want 
to display the asymptotic behavior of solutions at very large values of mi/2 and mo, but 
will do so in Section ^ where we study the implications of this work for SUSY searches at 
accelerators. We will see that, for some choices of rrit, tan f3, and Aq, both mo and mx/2 are 
bounded from above by purely physical criteria, and no fine-tuning constraint is needed. 

8.2 constrained Minimal PArameterS Space (COMPASS) 
8.2.1 General Properties 

We now focus on the region of the parameter space consistent with all the constraints listed 
above. This region certainly meets our expectations for where SUSY might be realized 
because the gauge couplings unify there, correct EWSB takes place, and the experimental 
and cosmological constraints are satisfied. In this constrained region of parameters (COM- 
PASS) we now analyze the various relations that result between the SUSY spectra and 
the implications for SUSY searches. Next, we will study what additional restrictions are 
implied by imposing the dark matter constraint. 
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Mass Limits 
(GeV) 


COMPASS 


CDM 


MDM 


lower 


upper 


lower 


upper 


lower 


upper 


h 


61 


79 


62 


73 


61 


71 


A 


635 


1934 


691 


1340 


658 


1115 


gl 


183 


595 


241 


403 


208 


332 


eR 


111 


408 


190 


267 


141 


207 


n 


110 


407 


190 


267 


140 


207 


T2 


183 


595 


241 


403 


208 


332 


VL 


176 


592 


236 


400 


202 


328 


UL 


550 


1621 


571 


1129 


559 


943 


UR 


530 


1549 


552 


1082 


539 


905 


h 


342 


1199 


354 


810 


347 


660 




607 


1546 


620 


1112 


612 


948 


=LSP 


97 


356 


97 


233 


97 


189 




182 


669 


183 


440 


182 


356 


xt 


180 


668 


182 


440 


180 


355 


9 


596 


1780 


597 


1234 


598 


1028 



Table 4: The lower and upper limits for the case mt = 145 GeV, tan/3 = 1.5, Aq/itio = 0, 
and sgn//o = ~1 (Fig- 0) ^'^^ solutions in COMPASS (with no fine-tuning cut 

satisfying / < 50 imposed), and for the subset of solutions selected by either the MDM 
or CDM constraint. Because of the finite-size grid in our numerical sampling the limits 
presented here could be somewhat relaxed and should be treated only as indicative. 

Several typical examples of solutions resulting from our analysis are presented in more 
detail in Figs. |7[]9| and in Tables In the graphs we show the typical ranges of several 
interesting parameters. In the tables we display the lowest and largest values of various 
masses selected after scanning all the choices of parameters compatible with COMPASS. 
(The ranges selected by DM, presented in the last two columns, will be discussed shortly.) 
We see that the allowed mass ranges are rather broad and typically allow for masses as light 
as, or not much heavier than present experimental limits. 

On the other hand, we see that, without constraining mi/2 and ttiq from above by the 
fine-tuning constraint, all the masses can (for some rrit and tan /3) become very large, with 
the squark, gluino, and heavy Higgs bosons {H, A^ and H^) typically being the heaviest 
and the sleptons, charginos, and neutralinos being significantly lighter, except for mo large 
and mi/2 ^ 0{mz) where mj k, m~ ^ m~. Very large values of rriQ ^ ^1/2 ^.nd large 
values of mi^2 ^ ^0 are typically disallowed by ^x^o < 1' charged LSP, color breaking 
(tachyonic ti), and no EWSB. But for many choices of parameters, one can only exclude 
both large mi/2 aiid mg by imposing the fine-tuning constraint. Thus we see that without 
further constraints or criteria, COMPASS still allows for a wide range of SUSY masses, 
though these masses are correlated in very specific ways. 

An important quantity in the MSSM is the Higgs/higgsino mass parameter ^u. In con- 
trast with "^0) and A* does not break SUSY and therefore a priori it could be 
much larger than mz- Similarly, while supergravity suggests a value of order mz for mo and 
"^1/2) if generically does not say anything about the origin of /xo = /x(Mx). On the other 
hand, phenomenologically, it would be very surprising if one of the defining parameters of 
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Mass Limits 
(GeV) 


COMPASS 


CDM 


MDM 


lower 


upper 


lower 


upper 


lower 


upper 


h 


91 


113 


96 


112 


94 


108 


A 


209 


1773 


346 


1402 


314 


1115 




118 


1208 


202 


1047 


175 


1022 


eR 


84 


1070 


162 


1015 


127 


1007 


Tl 


81 


1067 


160 


1012 


125 


1004 


T2 


120 


1206 


203 


1046 


176 


1021 


l^L 


89 


1205 


187 


1044 


157 


1019 


UL 


326 


2175 


480 


1965 


417 


1368 


Ur 


318 


2095 


465 


1877 


407 


1310 


h 


196 


1615 


310 


1552 


250 


1064 




408 


2011 


538 


1877 


481 


1333 


=LSP 


29 


437 


76 


435 


59 


286 


xi 


63 


801 


122 


789 


106 


521 


xf 


51 


800 


112 


788 


101 


520 


9 


294 


2153 


502 


2151 


419 


1491 



Table 5: The same as in Table ^ but for mj = 145 GeV, tan/3 = 5, Ao/mo = —1, and 
sgn//o = -1 (Fig. D- 



Mass Limits 
(GeV) 


COMPASS 


CDM 


MDM 


lower 


upper 


lower 


upper 


lower 


upper 


h 


113 


131 


116 


125 


114 


119 


A 


532 


1502 


564 


1020 


532 


828 




244 


1069 


244 


1011 


244 


832 


eR 


167 


1023 


167 


1004 


167 


824 


Tl 


144 


980 


144 


960 


144 


788 


T2 


250 


1051 


250 


991 


250 


816 


VL 


230 


1066 


230 


1008 


230 


828 


UL 


641 


1681 


677 


1156 


641 


931 


Ur 


631 


1611 


654 


1110 


631 


924 


h 


441 


1302 


501 


883 


464 


607 




584 


1579 


687 


1117 


605 


814 


Xi =LSP 


28 


353 


34 


232 


34 


152 




51 


657 


62 


432 


62 


281 




50 


657 


61 


432 


61 


281 


9 


207 


1812 


249 


1257 


249 


874 



Table 6: The same as in Table | but for mt = 170 GeV, tan^ = 20, Ao/mo = 0, and 
sgn//o = -1 (Fig. I). 
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the MSSM were much larger than others. In our analysis fx is determined by the other 
input parameters and the adopted constraints. We find |/i| broadly in the range of values 
spanned by either or tjiq. Two typical patterns can be identified. In the cases when the 
constraint from EWSB does not exclude the upper left-hand part of the {mif2,'mo) plane, 
we find 1^1 ~ mi/2 for ^ ^0 and |//| ~ mo for niQ » mi/2- Otherwise, \fi\ ~ mi/2 

for small thq but slowly decreases as mo grows. Overall, the values of resulting from the 
analysis are closely related to mi/2 and mo and only for such values of fi does the CMSSM 
appear to be self-consistent. 

One important consequence is that the lightest neutralino x is in most cases gaugino-like 
(more specifically, bino-like) ||2^. (Figs. and Tables ^jG] show typical neutralino mass 
ranges and compositions; see also Fig. ^2|). This is quite a remarkable theoretical prediction 
of the CMSSM in light of the fact that a bino-like neutralino has been selected theoretically 



as the unique attractive candidate for (neutralino) dark matter |84, 91]. Notice also that, 
while X is typically at least 80% (and in most cases 90%) bino, it is never a pure bino state. 
Various analytic approximations for ^x^o and related bounds derived for a pure bino may 



thus be misleading [96|. 



It is worth noting that the LSP has typically a dominant bino component because \fj,\ 
almost always comes out somewhat larger than M2 ~ 0.8mi/2 (compare Fig. [3l| ). The 
composition and scaling properties of the neutralinos and charginos in the plane (/U,M2) 



have been well-understood [91, 92, In particular, for \fi\ > M2 the lightest neutralino is 
mostly gaugino-like (in fact, bino-like; see, e.g., Fig. 1 and the discussion of gaugino purity 
in Ref. [^]). For gaugino-like LSP the masses roughly satisfy the relations 

m^ ~ Ml ~ O.5M2, (32) 

m^o ~ m^± ~ 2m^, (33) 

(These approximations improve as ^ M2.) It is important to note that in this approach, 
these relations are characteristic to most solutions in COMPASS, and do not come from 
GUT-dependent constraints, such as proton decay. 

In some regions of the (miy2)''^o) plane we do find LSPs with significant higgsino com- 
ponents. This happens for both mi/2 and mo small ( < 100 GeV), the region typically 
excluded by experiment. It also happens in relatively small regions close to where EWSB 
cannot be achieved. There |/u| is smaller than mx/2- In a few other cases we also find 
higgsino-like LSPs for larger mi/2 and mo. This happens for small rrit, tan/? well above 
one, and very large ^0 (e.ff-, for nit = 120 GeV, 3 < tan/3 < 20, Ao/niQ = 3, sgn/^o = ±1) 
in a relatively limited region of large mi/2 ~ n^o ^ 400 GeV disfavored by fine-tuning (com- 
pare Fig. [T^). Higgsino-like LSPs have been shown, however, to provide very little relic 
abundance [g^. For > mz,mw ,mt the x pair- annihilation into those respective final 
states {ZZ, WW, tt) is very strong [96|. Both below and above those thresholds, there are 
additional co-annihilation |^6| processes of the LSP with xf and which in this case are 
almost mass-degenerate with the LSP. Co-annihilation reduces ^x^o b^low any interesting 



level [|24| , |97| . Higgsino-like LSPs thus do not solve the DM problem. Except for those 



relatively rare cases, we find an LSP of at least 80% bino purity. 
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It is also interesting to explore what values of ^ at the GUT scale (/ig) result from the 
analysis. We choose to display it in terms of the ratio |^o|/?t^o- A very simple relation 
emerges: |/Uo|/?7io decreases from a few in the large and small thq region down to one 
or less in the opposite extreme. If, for some choices of parameters, the ratio falls down to 
zero, no proper EWSB occurs. 

One other parameter of the model is B, which does not run very much between its GUT 
value Bq and B{mz) (see Table ^). A typical tendency is for B/mo to grow with mi/2 and 
decrease with mo- We also show in Fig. ^ a scatter plot of Bq vs. Aq for mt = 145 GeV 
for all the solutions belonging to COMPASS. The value of Bq that we obtain as an output 
of our procedure rarely yields the relation Bq = Aq — tuq that is often imposed by other 
analyses. 

In this Section we have focussed mostly on the case mt = 145 GeV. Varying mt leads 
to significant modifications but the general features remain, as can be seen by comparing 
Figs. 0-^. We find for mt = 170 GeV that the constraint coming from imposing proper 
EWSB becomes much weaker. Similarly, the regions where the lightest neutralino is not 
the LSP become pushed towards even larger 772-1/2 • For a given point in the ("2-1/2) ^0) plane 
Osi'iTT'z) grows with mt mostly due to the sin^ dependence on mt (Eq. (|2^)), as discussed 
in Section |6.2| . 



8.2.2 Regions Favored by the Dark Matter Constraint 

We now point out the sub-region of COMPASS which is favored by the hypothesis that 
the LSP is the dominant component of either cold or mixed dark matter. As we discussed 



in Section [7.4.41 , there is now abundant evidence for the existence of DM in the Universe. 
The neutralino has become one of the most attractive candidates for DM. In the pure 
CDM scenario one expects the LSP relic abundance to be in the range given approximately 
by ([29|), while in the currently more favored mixed (CDM+HDM) scenario it should roughly 
satisfy the range (^0|) . 

Applying either ( |29| ) or (^) to the parameter space under consideration results in se- 
lecting only relatively narrow bands in the plane (7ni/2, ^n-o) whose shape and location vary 
with other parameters but typically correspond to both mi/2 and mo in the range of a few 
hundred GeV. (See Figs. 0-^.) Of course, they fall into the region constrained by the age 
of the Universe {^x^o < 

More importantly, requiring enough DM (i.e., taking lower limits in either ( p9| ) or (pC|)) 
typically leads to /oiiier limits on both r?T-i/2 and mo and, as a result, also on the SUSY mass 
spectra which are higher than in COMPASS alone. It is interesting that the mass ranges 



consistent with either (2£) or (|30| ) are typically less accessible at LEP II and FNAL. This 
can be seen by comparing the lower limits allowed by COMPASS with those selected by the 
CDM or MDM scenarios in Tables ^[-^. (See also Table § and the discussion in Section |9.9| .) 
For example, in the case presented in Fig. ^ and Table ^ applying the DM constraints causes 
the chargino xf and the sleptons to be completely inaccessible to LEP II and the gluino 
to be above the reach of FNAL. It also makes it harder to discover h and other particles. 
On the other hand, the DM constraint severely lowers the upper ranges of masses for all 
the particles making them much more likely to be accessible at future accelerators like the 
NLC or LHC. Prospects of searches for various particles will be discussed in more detail in 
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Section ^, and in particular the detectability of the hghtest Higgs as a function of LEP II 
beam energy will be analyzed in Section 9.2.2| . Here we only note that, with large enough 



y/s, h has a very good chance of being discovered at LEP II. 

While one might argue that the constraints (]29| ) or (^) do not carry the same weight as 
some other constraints listed above, they do reflect our current cosmological expectations 
and serve as a strong guide to those regions of the parameter space in which SUSY solves 
the DM problem. 

8.3 Effect of the Full Effective Higgs Potential 

The results from our analysis have also served to reinforce the need for using the full 1-loop 
effective potential in the minimization procedure |51]. We already argued in Section || that 



the 1-loop contributions to 14iiggs were important in order to stabilize the scale-dependence 
of the potential, but one can also see the net effect of using the full 1-loop Vniggs in our model- 
building results. As well, one can see the smaller role played by the non-leading contributions 
to Twiggs, that is, contributions not coming from the t - 1 splitting |5|, |5|] . In Fig. |l|, 

we have shown two plots of the 5 "^o) plane for the choice of input parameters as in 



Fig. ^. Fig. 15a shows the region of parameter space allowed after we have excluded the 
regions in which EWSB did not occur (labelled E), where the LSP was charged or colored 
(labelled L), and where the neutralino relic abundance would "overdose" the Universe 
(labelled A), for the renormalization group- improved tree level potential only. On the other 
hand. Fig. |l5|b shows the parameter space available for the same choices of parameters, but 
now with the renormalization group-improved 1-loop effective potential with leading terms 
only. Notice that the regions in which EWSB did not occur have enlarged, taking over 
some of the regions which were excluded before on the basis of their LSP being electrically 
charged. However, the strong bound placed on the parameter space by DM constraints has 
considerably weakened, leaving the region in which mo — mi 12 large available, pending a 
fine-tuning cut. We note also that the effect of including the 1-loop contributions to Vniggs 
is negligible in the region mi/2 *C tuq favored by the proton decay constraint |19, ^] in 
the minimal SU{5) GUT. Finally, including all (leading and non-leading) terms does not 



modify the situation sizably as can be seen by comparing Figs. 15b and pa. The qualitative 



difference is extremely small, which we found to be a general result. 



9 Applications 

9.1 Overview 

The analysis described in previous sections has led us to a restricted parameter space for 
rrit, tan/3, mi/2, rriQ, Aq, and sgn/^o = ±1 in the CMSSM which we call COMPASS (see 



Section |8.2[) . Most previous studies of SUSY predictions have preferred to fix some of these 
parameters by assumptions and vary one or two, either with or without constraints. This 
is useful and interesting and can lead to instructive predictions, but there is always doubt 
about their generality. 

We have taken the alternative approach of studying the fully constrained parameter 
space described in Section ^. We know that any point in COMPASS is already guaranteed 
to have gauge coupling unification, a Higgs mechanism, all phenomenological constraints 
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satisfied, etc. We can then ask a variety of questions about the regularities of the resulting 
solutions, whether they have predictions of interest, and so forth. For example, we can ask: 
what fraction of solutions gives a spectrum of sparticles that can be detected at LEP II and 
FNAL (or any other present or future facility), and in what channels do we most expect 
to find sparticles? What do the solutions predict for BR(5 57), T{Z ^ bh)l What is 
il^/ig for the solutions? When new experimental or theoretical information is available it 
can be easily added to constrain the parameters further. In the following we describe a 
number of such results. More specifically, in this section we examine the solutions that pass 
all the theoretical, experimental, and cosmological constraints listed in Sections |7|-^, i.e., 
solutions in COMPASS. We impose two additional cuts. We keep only those solutions which 
require no large fine-tuning of parameters. We take / < 50 which roughly corresponds to 
the heaviest squark, gluino, and Higgs masses falling below 1 TeV, except for very large 
tan/? where m~ and m~ can be larger (see Section ^.1.8 ). We also impose the lower bound 



BR(6 57) > 1.5 X 10~^ |70|] (see Section |7.2D . The solutions in this restricted set will be 
called "acceptable." Most of the results presented in this section have been derived with 
rrit = 145 GeV, but in some cases we consider other values of mj. 

Recall that our constraints do not require a detailed knowledge of the physics at the high 
scale. Our parameter space is intended to be the most general one which is independent of 
multifarious GUT scenarios. It is for this reason that we do not impose a constraint on the 
lifetime of the proton. The proton decay constraints have been included first by Arnowitt 
and Nath [19|, and also by Lopez, et al. [30| mainly in an SU{b) GUT. They find a longer 
proton lifetime for smaller tan/3 and smaller mi so this region of the parameter space is 
enhanced for them. We will study the implications of adding assumptions about unification 
and a GUT group in the near future. 

We have also assumed a common scalar mass mo and a common gaugino mass mi/2, 
both of which can be relaxed, which we will consider in the near future. Keeping these 
comments in mind, we now discuss some CMSSM (Constrained MSSM) results. 



9.2 Higgs Physics 
9.2.1 What is m/j due to? 

In a supersymmetric theory with electroweak symmetry breaking the Higgs boson mass is 
calculable, and it is very interesting to ask what parameters in th theory play a role in 
determining the value of rrih. The tree level mass matrix for the two CP-even scalar bosons 
is 

2 ^ ( -B^i [3 + \{gl + gl)v'^ cos^ (3 Bfx - ^{gf + g^)v'^ sin (3 cos P \ , . 
y Bfi-l{gl+g^y sin (3 cos (3 - B cot (3 + ^gl + g^y sin^ (3 J ^ ' 

where = ^{gf + g2)v'^ and v'^ = v"^ + v^. If = or if w'^ = then this matrix has a 
zero eigenvalue. Thus in supersymmetry one cannot think of ruh as coming only from the 
Higgs self-interaction. Any interpretation is complicated since Bfi, tan /?, and v'^ are all 
involved. Furthermore, the one-loop effective potential can yield m/j significantly above the 
tree- level result p3[ ]. 

Haber has emphasized that there is a lower limit on nih', even if the tree level value 
is zero the one-loop potential generates a mass. He finds a lower value above 60 GeV, but 
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that assumes 1 TeV squark masses. We agree that there is a lower limit, but it is sensitive 



to squark masses, as shown in Fig. |16| where we plot the lower limit of ruh versus .yrnf^rnf^ 
which contributes the largest radiative correction to mh- The lower limit basically arises 
because of the way the EW breaking comes about in SUSY. Eq. ([l5| ) leads to a lower limit 
on -B//(= 7713) and thus a lower limit on m/i. 

Effectively, rrih arises from three sources: the product of SUSY parameters Bfj,, the 
value of the vevs (whose sum in quadrature is fixed numerically by mz), and the one-loop 
radiative corrections. To demonstrate how these sources of Higgs mass interplay, we show 
in Figs. 17 and |l^ plots of m/j vs. tan /3 and vs. ^/\B^\, where ruh is the full radiatively 
corrected Higgs mass. The upper limit on ruh is due to the usual argument that in SUSY 
the Higgs self-coupling is fixed by the gauge couplings with an additional contribution from 
the radiative corrections [p9|]. 

One can see a strong correlation between tan/3 and the allowed ruh- This is expected 
since the tree level upper bound for goes like |cos2/3|. In fact, we find from our solutions 
that, for mt = 145 GeV, tan /3 < 5 if m/j < 85 GeV. Therefore, if LEP H finds the Higgs 
then tan /? is constrained to be less than 5 for all our surviving solutions. Solutions with 
tan /? > 5 and rrih < 85 GeV are excluded mainly by one of three effects: (1) the chargino or 
sneutrino mass is too low; (2) the LSP is not the neutralino; or, (3) electroweak symmetry 
breaking does not occur. Fig. ^ shows the distribution of ruh for all acceptable models 
with mt = 145 GeV. 

9.2.2 Detection of the Higgs Boson 

Interestingly, we find essentially no solutions for which the h + A mode is detectable at 
^/s < 210 GeV (fewer than 0.1% of the solutions), since ruA is too large for all CMSSM. 

Almost all solutions have sin^(/3 — a) > 0.98, so the Zh cross section (which is propor- 
tional to sin^(/3 — a)) is not suppressed |^T|. Thus the experimental limit on the SM Higgs 
boson effectively applies to the h of the MSSM in all acceptable solutions. The current LEP 
bound is ruh ^ 62 GeV That sin^(/? — a) 1 and that tua is not small enough for the 
h + A channel to be accessible at LEP are related. 

A similar result holds for the tth coupling. It has a factor cosa/sin/3 which is within 
a few percent of one over essentially the entire set of solutions, so that methods to detect 
h by radiation off a top quark will work essentially as well for the SUSY h as for the SM 
Higgs. 



In Fig. 20 we show the percent of solutions with mt = 145 GeV for which h is detectable 
at a given ^/s. One can see that about 30% of the solutions are detectable when LEP 
energy increases to 178 GeV, increasing to about 75% if ^/s is increased up to 210 GeV; for 
the MSSM 100% is reached at ^ ~ 220 GeV. This hmit is weh known since in the MSSM 
the upper limit on m^ is about 125 GeV for mt < 150 GeV. 



A number of groups |100| have examined the detectability of at least one SUSY Higgs 



at LEP or SSC/LHC. They concluded that much of the complete parameter space could be 



covered, but not all. Results were often presented on a tan (3 vs. plot. In Fig. 21 we show 
where our constrained solutions appear on such a plot. We also mark the approximate region 
within which the detection of at least one SUSY Higgs was found unlikely |100| ] . Amusingly, 
about 2/3 of the solutions that do fall in the region are detectable at LEP II or FNAL in 
some other channel. 
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We also present in Fig. 22 a scatter plot of rrih vs. rriA for all acceptable solutions for 
rrit = 145 GeV. The distinct branches seen in the graphs correspond to different choices of 
tan /5. One can see how nih grows with rriA and tan /3. The dependence is smeared to some 
extent by varying all the other parameters of the model. 

9.3 What Is the Origin of m^? 

We have remarked on several aspects of the role of mt in other sections. Here we discuss 
briefly the questions of what contributes to the mass of mt- Of course, the main question 
is why nit ^ "T-fo and how that is answered in supersymmetry. More explicitly, once we are 
below the scale where EW breaking has occurred, one can write at a scale Q, 



mt{Q) = ht{Q)v{Q) tan P{Q)/^l + tan^P{Q), (36) 
for the top quark running mass, and likewise for the bottom, 



mt{Q) = h{Q)v{Q)/,Jl + tan2/3(Q). (37) 

If the SU{2) symmetry were not broken here, we might expect mt to be only a little larger 
than mi). To understand the effects that can enter, we can take the ratio at mt'- 

miQ = mt) htijnt) 
mb[Q = mt) hb{mt) 

and we can define 

_ ht{mt)/hb{mt) 

— T. 77 ^ ' 

where ht^bo = ht^biMx)- Then finally 

^^ = r^^tan/3. (40) 
mb{mt) hbo 

Thus the large ratio mt{mt) / mb{mt) ~ 50 could be due to any of three factors: tan/3, the 
ratio of the Yukawa couplings at the high scale, and/or the ROE running of the Yukawa 
couplings, as expressed by the value of r. 

The RGEs of the top and bottom Yukawa couplings are 

dht _ ht f 8 2 3 2 13 2 , ^,2 , 1^2 



dlnQ 87r2 V 3 



9l - lal - + 3/.? + \hi) , (41) 



dhb hb f 8 2 3 2 7 2 , o,,2 1,,2 li,2 



o . -gi- -59 9^ + 3/ih + -hf + -K] . (42) 

dlnQ 87r2 V 3^^ 2^^ 30^^ ^ 2^ 2^J ^ ^ 

Since the gf and /i^ contributions are numerically very small, we see that if hb ~ ht at the 
high scale, they will run down together, in which case the large value of mt is generated 
predominantly by large Vu, and necessarily tan/? ~ mt/mb is large. For ht larger than hb 
at the high scale, mt increases relative to mb from the running, and the physical mt is 
reached with a smaller tan /3. To illustrate these effects we show in Fig. |2^ a scatter plot 
of hto/hbo vs. tan/3 for the constrained solutions, all with mt = 145 GeV. We see that 
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{hto/hiio) tan/? ~ 60 is a good approximation to the results except for very large and small 
tan /3, so that the large value of rrit cannot be interpreted as coming from the running; it 
must be input, either as a large GUT-scale ratio /ijo/^w oi' &s a large tan/3. In particular, 
if h is discovered at LEP178, then the large mt must be due to the top Yukawa at the GUT 



scale (compare Fig. 17). 



9.4 BR(b^s7) 

The recently reported upper bound BR(6 57) < 5.4 x 10~^ (see Section |7.2D has spurred 



an increased interest in predictions for 6 ^ 57 in SUSY. Barbieri and Giudice [73| have 
reminded us that in the limit of unbroken supersymmetry the MSSM prediction (including 
the Standard Model part) is zero due to a theorem of Ferrara and Remiddi [|101| , and 
they have shown that SUSY solutions will give reasonable values for this rate. Garisto 
and Ng ||76|, and others, have also done a general SUSY analysis of the implications of 
BR(6 ^ S7). 

In Fig. ^ we show a histogram of BR(6 57) for all the solutions. We have checked 
that the solutions in the peak come from all over the parameter space and in no sense 
represents a decoupling region. For a typical solution the magnitudes of the W-t loop, the 
H^-t loop, and the x^-'t loop contributions are all about the same, with the W-t and H^-t 
loops having the same sign and the x^-t loop having the opposite sign. We see that the 
CMSSM naturally produces solutions in the right range. These results are predictions in 
the sense that here we impose no constraint on the model space from 6 — > 57 data. (For 
other uses of the model space we cut at the upper and lower limits indicated in the figure, 
so that our model space does include the BR(6 57) constraint in general, except for the 
present discussion.) 

Some authors 173] have claimed that this decay strongly constrains charged Higgs boson 



masses. To show that there is no strong constraint in the CMSSM, we show in Fig. 25 a 
plot of BR(6 — > 57) vs. mjj+ for the solutions in the region between the upper and lower 
limits. That is, every mjj+ in Fig. ^ gives a BR(6 57) consistent with experiment. 



Finally, it is interesting to look at BR(6 — > 57) vs. tan/3 in Fig. 26. Smaller tan/? values 
concentrate somewhat in the allowed region, though acceptable solutions occur at any tan (5. 
In the large tan /? region the chargino contribution can be quite large and negative |^, |7^ . 

9.5 Detection of SUSY at LEP II and FNAL 

There are a number of possible ways to detect supersymmetric partners at FNAL or LEP II. 



We estimate that about 32% of all CMSSM acceptable solutions (see Sec. |9.lD have either a 
superpartner or the light Higgs boson detectable at LEP with y/s = 178 GeV and 500 pb~^, 
or at FNAL (with, say, 500 pb~^ integrated luminosity), or both {h will only be detectable at 
LEP II, not FNAL). We include in this sample all solutions with mt = 145 GeV, mo < 1 TeV, 
^■1/2 ^ ITeV, and / < 50; that is conservative, giving q and g masses over 2 TeV for the 
largest mo, mi/2- The solutions with large mo, mi/2 are generally not accessible at FNAL 
or LEP II, but they also do not increase the number of acceptable solutions rapidly enough 
to dilute the 32% result even if larger mg, mi/2 were to be included. The results do not 
vary rapidly with mt- These numbers are for solutions with ^x^o ^ 



38 



At LEP II 24% of the acceptable solutions allow detection via one-sided events, e+e^ 
X2X1) where Xi = X) followed by X2 ~^ ^~^^~Xi- And 18% will have pair production of 
the lightest chargino, and 8.6% detection of h. There is overlap, of course, and 30% of all 
solutions are detectable at LEP II. Selectrons are detectable in 3.6%, a light stop 0.7%, and 
h + A in 0.06% of the solutions. 

An interesting way — perhaps the only way at LEP — to determine if /i is a SUSY Higgs 



is to measure the cross section for e"'"e~ — > /i + nothing |102| . In the Standard Model this 
entire cross section should be from e^e^ h{-^ bb) + Z(-^ vU) and will be very accurately 
known. In SUSY there is also a contribution from e^e~ — > ^{-^ ^+Xi) + Xi) unfortunately 
this contribution is larger than 10% of the SM cross section in only 0.6% of the solutions 
at LEP 178 but would be the "proof" of SUSY in these cases; at larger y/s the fraction of 
solutions where this effect could be observed increases rapidly. 

As an illustration, we show in Fig. ^ what fractions of the ("10)^1/2) plane would be 
constrained by SUSY searches at LEP II if nrit = 145 GeV. The kinematic criteria that we 
use to determine the detectability at LEP 178 are m^± < 85 GeV, mj < 85 GeV (for any 
sfermion), m^o + m^o < 170 GeV, rrih + rriA < 170 GeV, and nih + mz < 170 GeV. We 
also require that any event-signature lepton have energy above 5 GeV or any quark have 
energy above 10 GeV. The regions marked by crosses (empty boxes) will always (never) be 
accessible to LEP II for any combination of input parameters. Filled boxes mark the regions 
accessible for some combinations of parameters. In window (a) we show the combination of 
possible SUSY searches at LEP II by applying the criteria listed above. In window (b) we 
show the same for the chargino Xi^ alone, and in window (c) for the lightest Higgs assuming 
rufi < 80 GeV. {h is a very SM-like Higgs.) Finally, window (d) shows how much larger 
a region would be explored by searching for h up to 110 GeV. Remember that very large 
values of mi/2 a-^e disfavored by the fine-tuning constraint (compare, e.g., Fig. 13). 

At FNAL gluino detection will occur in 11% of all solutions, squark detection in 5%, 
detection of Xixt in 25%, xfxi in 14%, X2X1 in 24%, X2xf in 12%, 1 1 + 1 1 in 4%. These 
combine to make 26% of all solutions being detectable at FNAL. The FNAL/LEP overlap 
is large, so combining them only increases the percentage of solutions detectable at FNAL 
or LEP II to the above mentioned 32%. 

The kinematic criteria that we use to determine the detectability at FNAL are mff± < 

mt-5GeV, m~ < 300 GeV, m + < 85 GeV, m~< 300GeV, andm^o+m^o < 170 GeV. We 
y Xi y Al A2 

also require that any event-signature lepton or quark have energy above 15 GeV. Some of the 
percentages listed above for the detectability of different channels will decrease, particularly 
at FNAL, when detection efficiencies and cuts to reduce background are considered. But 
with sufficient luminosity and sufficiently good detectors the above numbers should be 
approached. We are presently undertaking full simulations of signals and backgrounds to 
determine reliable signatures and strategies. There have been previous studies of FNAL 
and LEP II detectability in some depth [|103(| . Our only advance so far over some of these 
is that their conclusions are based on a parameter space some parts of which are excluded 
because the various constraints are not satisfied. 

We note that at FNAL, for gluinos lighter than 300 GeV, 80% of solutions have squarks 
heavier than gluinos, so that the appropriate way to simulate g detection is to take m~ > m~ 
in the first approximation. 

We have also investigated our CMSSM parameter space to see how the top quark search 
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at FNAL could be affected by supersymmetry. We find that for mt = 170 GeV approx- 
imately 3% of all acceptable solutions kinematically allow one or more of the following: 
t bH^ decay, t tiXi decay, or g tit with g < 250 GeV. (For smaller mt this fraction 
is reduced.) Any one of these kinematic possibilities can significantly alter the kinematic 
analysis and/or the effective rates (after cuts) of top quark production. Once the top physics 
at FNAL settles into place if none of these is observed then the parameter space is reduced 
a few percent. 

Some reduction in detectable solutions would occur if the xf ~ Xi mass difference were 



so small that the resulting lepton or jet from Xi decay were too soft to detect. Fig. 28 



shows a plot of m ± — m o vs. m ± from which we see that most solutions have no problem 

Xl -^1 Xl 

here; and Fig. |2^ shows the energy of the lepton or jet from xf decay. 

We will report a study on how effective higher energy linear colliders will be at studying 
SUSY for the constrained model space later. For now we note that NLC with ^/s = 350 GeV 
will be able to detect h and at least one superpartner for about 75% of the constrained 
solutions; this number grows to about 97% as ^/s grows to 500 GeV. 

If h is not detected at LEP II once ^/s ~ 210 GeV, most but far from all solutions will 



be excluded, particularly if mt ~ 145 GeV. Fig. 19 shows a histogram of m^ values and one 
can see that most solutions are below 110 GeV. In solutions with a Higgs sector extended 
beyond that of the minimal one there is still an upper limit on mh, but it can be as large 



as about 146 GeV M 



Overall, we conclude that, while not finding superpartners at LEP II and FNAL does 
eliminate nearly a third of the parameter space, it will still leave many possibilities open. 

9.6 What If FNAL and LEP II Do Not Detect a Sparticle? 

If superpartners are not detected at LEP II or FNAL then much of the low (mo, ^1-1/2) 
parameter space can be excluded. Bounds on mi/2 are determined mainly by the bounds 
on the gluino. The gluino mass is related to mi/2 by 

mi/2 = (43) 

We find the lower limit of ^~ = ax/as(m,~) to be (compare text below Eq. (p2|) ) 

e~ > 0.36. (44) 

So if m~ is determined from experiment at FNAL to be greater than 300 GeV then 

mi/2 > (0.36) (300 GeV) = 107 GeV. (45) 

Bounds can also be placed on mx/2 by direct searches on the lightest chargino. By taking 
the square (MM-^) of the chargino mass matrix, we find that 

< i^±m\i2 + 2mwCos^/3. (46) 



Xl 

where 



^2 16 
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So if m ± > Mw then the tan /^-dependent bound for rrii/o is 



"1-1/2 > "^VK^Y 1 — 2cos2/3 



(48) 



For high tan /? this bound (~ 100 GeV) would become comparable to the resulting bound 
on mi/2 from the gluino search. In any case, bounds on mi/2 are obtained straightforwardly 
from direct sparticle searches. 

It is more difficult to put bounds on ttiq. Although the squarks and sleptons gain 
mass from uiq they also have contributions from the gaugino masses (compare Eq. (p^)). 
Therefore bounds on mo from direct searches of squarks and sleptons must be presented 
as a function of mi/2- Letting / represent any squark or slepton and ifij the lower mass 

bound of / we can write the inequality relation that constrains the {mQ,mi/2) plane: 



2 I u 2 ^ ~ 2 

™0 + ^7^,«"^l/2 > 



For example, for / = e^. 



m'j =p cos 2/3 



9 u 1 
— <b~ < -. 

20 2 



rpfL,R 



QfL,R 0^ 



(49) 



(50) 



The upper and lower bounds on are good to within 5% accuracy. As can be clearly seen 
from Eq. ( ^9[ ) the slepton bound translates to an ellipse (first quadrant) in the (mo, 1^1/2) 
plane where all (m-O) '^1/2) inside the ellipse are excluded. 



9.7 What If FNAL or LEP II Does Find a Sparticle? 

If FNAL or LEP II discover one or more sparticles we would like to extract from this 
the GUT-scale Lagrangian. That is, we would like to extract the super symmetric input 
parameters mo, rni/2, Aq, tan/3 and sgn^o from all the observables that are sensitive to 
them. (We assume that we know m^). In this section, we briefly discuss how this should be 
done. It turns out that this is more difficult than it first appears; analytic methods are of 
limited applicability. 

The equations used in the last section to put bounds on mo and mi/2 given bounds on 
the gluino and sfermions can also be used to pin down mo and mi/2- In general, each fixed 
value of an observable, whether it be a mass or a cross section or an asymmetry or any- 
thing else, generates a hypersurface in input parameter space (mo, mi/2) tan/3, ^O) sgn/Uo). 
(The gluino mass, for example, would predominantly determine mi/2 np to additional un- 
certainties discussed below.) The effective dimensionality of this surface is determined by 
how many input parameters significantly affect the value of the observable. Determining 
this hypersurface is very difficult because of all the nonlinearities in relating low energy 
parameters to the input parameters through self-consistent solutions of many RGEs. In 



the last section and in Refs. [105| there are some simple equations that allow us to estimate 
functional relations between input parameters and sparticle masses. Even though these 
simple relations are often a good approximation to the full analysis we must keep in mind 
that every observable necessarily depends on all input parameters, though with varying im- 
portance. For example, a precise determination of mi/2 from m~ is limited because as{m~) 
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depends in general on all sparticle thresholds (at 1-loop) 



m~= —mi/2 = /(mo,mi/2, Ao,tan/3,sgn/io). (51) 

^ ax 

Likewise, precise gaugino contributions to sfermion masses require detailed knowledge of 
the sparticle thresholds which introduce all input parameters into the final determination of 
the sfermion masses. Methods based on analytic expressions are clearly limited. They may 
suffice to provide rough estimates of input parameters, especially soon after the discovery 
of a given particle, when the experimental bounds will be still large, but they cannot be 
improved to accommodate a complex analysis based on several well- measured observables. 

Our approach to the problem of determining input parameters from low-energy observ- 
ables does allow for such improvements. We explore wide ranges of the input parameter 
space and let the computer do the work. We have really already employed this technique 
to generate the COMPASS of the CMSSM. We have cut hypersurfaces in input parameter 
space with Higgs and sparticle mass bounds, BR(6 — > s-f) limits, ^x^o bounds, fine-tuning 
limit, as{mz), etc. However, highly precise measurements of observables especially sensitive 
to supersymmetric loop corrections, or direct measurements of sparticle production, will pin 
down the input parameters which determine those low energy observables. 

We present in Table |^ an example of how several experimental measurements can shrink 
the input parameters space to a "point" which would in turn allow us to predict all other 
observables (other sparticle masses, cross sections, etc.). The first row in Table ^ lists 
the range allowed for each input parameter given detection of just the lightest Higgs (m^ 
determined from experiment to be, say, ruh = 80 ± 5GeV). The value of fj, is evaluated 
at mz and its sign is the one of sgn/ig. We see that just knowing the mass of the lightest 
Higgs alone clearly does not significantly constrain any of the input parameters. The next 
row in Table |^ lists the range allowed for each input parameter given detection of just the 
lightest chargino {xt determined by experiment to be, say, m ± = 70 it 5GeV). Here 
again knowledge of just one mass, the lightest chargino, does not significantly constrain 
input parameter space. In the next three rows we list the range of input parameters given 
detection of just mr {mr = 110 ± lOGeV), just g {m~ = 210 it lOGeV), and just ei 
(mg^ = 95 lb 5GeV). The next row assumes detection of h and Xi ■ Notice how the 
combination of these two masses restricts mo, mi/2 and tan/3 far beyond what knowledge 
of each mass can do individually. 

As we progress down the rows of the table with each subsequent row assuming more 
and more detected particles, the input parameters become more and more constrained until 
mo, mi/2, tan/3, and sgn//o are determined precisely at the level of our numerical sampling. 
Only Aq remains stubbornly undetermined though it is better constrained. This is a general 
rule about Aq: few observables are very sensitive to it. Observables which depend on third 
generation sparticle left -right mixing are most sensitive to Aq. 

As this example shows, by generating many self-consistent solutions and "filling up" 
input parameter space with them, we have the means by which to use all observables 
simultaneously to constrain input parameter space. This approach of generating solutions 
with the most precision possible, calculating observables without untrustworthy simplifying 
assumptions, and then simultaneously comparing all the generated solutions with all the 
calculated observables is a powerful way to analyze and constrain minimal super symmetry. 
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Detected Particle(s) 


mo 


■mi/2 


Ao/mo 


tan/3 




h 


< 671 


< 549 


unbounded 


unbounded 


-472 - 253 




unbounded 


74 - 450 


unbounded 


< 40 


-477 - 332 


h 


< 549 


61 - 202 


unbounded 


< 20 


-428 - 420 


9 


unbounded 


74 


unbounded 


> 3 


-460 - 205 


EL 


< 74 


61 - 111 


unbounded 


< 5 


-237 - 172 


h, Xf 


< 368 


74 - 136 


unbounded 


< 3 


-414 - 169 


h, xf , ii 


< 302 


74 - 111 


-3.0 - -1.0 


< 3 


-414 - 140 


h, xt' 9 


< 136 


74 


-3.0 - -1.5 


3 


120 - 140 


h, Xi, ti, 9, eL 


61 


74 


-2.5 - -2.0 


3 


120 - 130 



Table 7: Table demonstrating how input parameter space (mo, ^O) tan (3, and sgn ^o) 

is constrained by detection of particles. The initial ranges of input parameters are listed in 
Sec. 3.1. For this example we assume detection of these particles to mean ruh = 80ib5GeV, 
m ± = 70 ± 5GeV, mr = 110 ± lOGeV, m~ = 210 ± lOGeV and mgr = 95 ± 5GeV. If 
more than one particle type is listed in the "Detected Particle(s)" column then the range for 
each of the input parameters is found from knowledge of all listed particle masses. Keep in 
mind that the ranges of parameters listed in the table are values obtained on our numerical 
sampling grid and therefore have errors associated with them corresponding to the grid 
spacing. For example, when we quote mi/2 = 74GeV we really mean that we find no 
acceptable solutions with — 61 GeV (the next lowest mi/2 value on our grid) and no 
acceptable solutions with ^ 91 GeV (the next highest mi^2 our grid). If our grid 
were very fine-grained, then we could quote ranges that would more accurately reflect how 
well the parameters were determined, and that reflect the experimental errors better. 
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It is this approach that will quickly enable us to add better measurements of observables 
and any new observables to the CMSSM, including possible announcements of sparticle 
detection. With our method we can go directly from data to the parameters of the effective 
Lagrangian at the unification scale. 

9.8 The r(Z hi) Partial Width 

Precision measurements at LEP II currently show a slight deviation from the Standard 
Model value for the extracted value of r^g/Fhad- The current LEP average [^] is 

Rb = r-^ = 0.2200 ± 0.0027. (52) 

Phad 

Several groups |106| ] have studied the loop corrections to this partial width in the Standard 



Model. The SM prediction for is heavily dependent upon the value of mt, but |46| the 
predicted value is 1.5a lower than the measured value for nit = 145 GeV and it gets even 
lower for higher mt- 



The super symmetric contributions to this decay width have been calculated in Refs. [107 



108(1 . We use the equations of Ref. |10q| to calculate the supersymmetric contributions to 



Rb within the CMSSM. We perform scalar integral reductions and numerical calculations 



from Ref. M 



In Fig. 30 we plot the histogram of all acceptable solutions with mt = 145 GeV and 
0.16 < ^x^o < 0.33. Notice that within the CMSSM the supersymmetric contributions 
tend to increase This increase in Rj, is mainly due to light ti and Xi^- Interestingly, if 
we require the calculated R^ to be within la of the measured value, then approximately 
75% of the resulting solutions will be detectable at LEP II and 83% of the solutions will 
be detectable at FNAL. The channel which is by far most detectable at LEP II for these 
solutions is XiXi production. At FNAL many possible channels for detection of sparticles 
are allowed: g, xfxi, xfxi, xfxt and XiX2- 

9.9 Neutrahno LSP as Dark Matter 

As we have already seen in Section |8|, the relic abundance ^x^o °f lightest neutralino is 
an important quantity which plays a significant role in constraining the parameter space of 
the CMSSM. If the Universe is at least 10 billion years old then ^x^o ^ (and the Universe's 
age of 15 billion years or more gives ^x^o ~ 0.25), see Section 7.4.3[ We have also seen 



in Section |8.2| that over significant regions of the model parameter space the LSP provides 
enough dark matter in either the CDM (0.25 ^ ^^x^o ~ 0-^; (H)) currently more fa- 
vored MDM (0.16 < Qx^o ~ 0-33, see (^) scenarios (see Sec. ffl^ ). This is a remarkable 
property of the CMSSM given the fact that, unlike in the case of many phenomenological 
quantities, calculating ^x^o involves elements of the physics of the early Universe and a pri- 
ori the resulting predictions for the LSP relic abundance could be completely incompatible 
with the expectation of low-energy SUSY. 

In this section we provide some more insights into the cosmological properties of the 



neutralino LSP. First, for the restricted set of "acceptable solutions" , as described in Sec. 9.1 



and for mt = 145 GeV we show in Fig. 31 a scatter plot in the plane (/i, M2). Notice a large 



concentration of solutions below the diagonals M2 = corresponding to the LSP being 
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mostly gaugino-like (see also the discussion above and below Eq. ^). This property is even 
more explicitly pronounced in Fig. 32 where we show for nit = 145 GeV and tan /3 = 10 a 
scatter plot of the LSP wavefunction (Zf^ for the bino, etc.) with an additional constraint 
0.16 < ^x^o < O-^"^ (MDM) imposed. Notice that the LSP is mainly a bino, as expected |]9l| . 
This has been already illustrated for a few specific cases in Figs. |^|^. A very similar plot 
can be made for any value of tan (3 and changes little for the CDM scenario. Without the 
MDM (or CDM) constraint imposed, we find in Fig. 32 also some points with somewhat 
smaller bino purity, with the largest concentration however still remaining at large Zf^. 

We also examine the predictions for the LSP relic density resulting from our restricted 
set of acceptable solutions defined in Sec. 9T, corresponding roughly to our expectations 
for low-energy SUSY. In Fig. 33 we plot a histogram of ^x^o otherwise acceptable 

solutions with mt = 145 GeV. Notice that there is a strong peak at ^x^o ~ ^-^ suggesting 
that in the CMSSM the MDM scenario is somewhat more favored relative to the CDM 
one. Within the framework of CMSSM we can view this result one of two ways, either as a 
prediction for 17^ (given /iq) for solutions which satisfy all other criteria, or that cold dark 
matter puts a severe constraint on the CMSSM if we demand that the LSP contributes 
most of the (cold) dark matter needed in the either the MDM or CDM scenarios. Both 
viewpoints are quite constraining: the first viewpoint for LSP (cold) dark matter and the 
second for the parameter space of the CMSSM. 

Finally, in Fig. 34 we plot ^x^o '^s- ^i/2: "^O) tan/3, and //. (In these graphs we 
lift the constraint ^x^o ^ Notice that the ranges of Qx^o favored by both MDM 
and CDM generally select both and \fi\ (which is an output parameter in our 

analysis) in a broad region of a few hundred GeV. For small 1711/2, rno ~ 100 GeV the relic 
density is too small because some sleptons are rather light there (roughly less than 100 GeV) 
which enhances the t-channel pair-annihilation xx ~^ ff 101 • Also notice that large tan /? 



produces more solutions with low ^x^o than do the solutions with low tan /3 which is due 
to a very strong enhancement in the LSP pair-annihilation cross section caused by the 
exchanged pseudoscalar A |Q. (The coupling Axx scales like tan/3.) 



9.10 The SUSY Spectrum 

By varying our input parameters over wide ranges of values, we can consider what ranges 



of sparticle and Higgs masses one should expect from the CMSSM. In Fig. 35 we present 



a scatter plot of mass vs. particle type for all the acceptable solutions in CMSSM with 
mt = 145 GeV. Little changes significantly with varying mt- In the next section, we will 
further prejudice our parameter space in order to find sample solutions that are preferred 
theoretically. 



9.11 Physics Prejudice Enhancement of Part of Model 
Space? 

In this section we apply some experimental and theoretical prejudices to the acceptable 
solutions. For example, solutions with BR(6 57) ~ 5.4 x 10~^ are less favored than 
solutions with BR(6 — > 57) « 3.5 x 10~^. Furthermore, theoretically we prefer solutions 
with lower fine-tuning. These are just two examples. Other prejudices that we can apply 
on the solutions are mb{Mx)/mT-{Mx) — 1, the LSP providing the right amount of cold 
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dark matter, and as{mz) close to its experimental central value. As can be gleaned from 
previous sections, some of these prejudices work against others in some respects. 

We have attempted to select a subset of all the solutions which are most likely to satisfy 
all (or most) of the above prejudices. We do that by effectively "squeezing" the solutions: 
select a preferred value for each constraint above and reducing the errors (or allowed region) 
by a factor of 2. In analyzing this squeezed set we find that the fraction of solutions which 
are detectable at LEP II and FNAL goes down by about a factor of 2. That is, only about 
18% of this set of solutions is detectable at FNAL or LEP II, whereas in the full set the 
fraction is 32%. We therefore find that the set of solutions which best satisfies our current 
experimental and theoretical prejudices are characteristically more difficult to detect than 
the full set of solutions allowed by current experiment. Table ^ presents three examples of 
such solutions. Solutions 1 and 3 are not detectable at LEP II or FNAL, but Solution 2 is 
in the chargino, h (LEP II) and gluino (FNAL) channels. We view this section as an initial 
attempt to add weighted physics criteria in order to select a part of the model space to use 
for other considerations such as phenomenological predictions or theoretical studies. 

10 Summary and Comments 

Encouraged by gauge coupling (grand) unification as implied by LEP, we have made an 
attempt to frame SUSY by re-considering Minimal Supersymmetry in the light of GUTs. 
We have parametrized the whole multitude of low-energy SUSY masses and couplings in 
terms of just five free parameters (and the sign of /i), including the mass of the top quark 
mt which will soon be known. We have demanded gauge coupling unification and proper 
electroweak symmetry breaking. In further reducing the allowed parameter space, we have 
included all the relevant experimental and cosmological constraints that can be imposed 
without choosing a specific GUT gauge group at the unification scale. We have not found 
the present experimental bounds on SUSY particle masses to be particularly constraining; 
in fact they are only beginning to limit the lower range of the SUSY masses. In contrast, 
rather robust cosmological constraints, like requiring that the Universe be at least 10 billion 
years old and that the LSP not be electrically charged, rule out large fractions of the 
SUSY parameter space. Furthermore, much more specific conclusions about the resulting 
SUSY mass spectra and properties can be drawn if one expects the neutralino LSP to be 
the dominant dark matter component (in either cold- or mixed-DM scenarios) in the fiat 
Universe. 

A number of groups have already reported studies along the same lines we follow. Our 
work is more comprehensive and complete in that more of the theoretical and phenomeno- 
logical constraints are included than in any previous work, and precision to the few percent 
level is required wherever appropriate in a fully consistent manner. We also include more 
applications than have been considered previously. 

It is only by combining all the constraints and exploring wide ranges of parameters that 
one is able to establish where SUSY might be realized. Remarkably, we find that SUSY is 
preferably realized in the range of Higgs, sfermion, and gaugino masses of several hundred 
GeV and below, with larger values sometimes allowed by our constraints but disfavored 
by too much fine-tuning. At this point one still cannot favor any range of mt, unless one 
insists on the rrib — m-j- unification which in most cases (but not always) implies a very heavy 
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Model Parameters 


Solution 1 


Solution 2 


Solution 3 


tau.i. 


10, 14.") 


1.5. 145 


."). 170 


mo, mi/2, Ao/mo 


247, 302, -2.5 


91, 111, 2.5 


111, 247, 2.5 


Bo/mo, B{mz) 


-0.67, -0.03 


2.84, 2.20 


0.02, -0.46 


Mo, IJ'imz), as{mz) 


394, 450, 0.124 


-214, -218, 0.127 


303, 304, 0.129 


ax, GcV 


0.041, 1.64 


0.042, 2.21 


0.041, 2.04 


h, H, A, 


116, 346, 345, 354 


62, 317, 305, 315 


113, 326, 324, 333 


CL, ML, TL 


328, 328, 324 


124, 124, 124 


211, 211, 211 


gr, JIr, tr 


276, 276, 268 


104, 104, 104 


152, 152, 152 




QIC Q1Q Q1Q 

oio, oio, oio 


11/1 11/1 11/1 

ii4, 114, 114 


1 07 1 07 1 (17 

ly ( , ly ly ( 


Ul, Cl, tL 


700, 700, 634 


283, 283, 299 


570, 570, 555 


Ur, Cr, tR 


677, 677, 620 


276, 276, 271 


551, 551, 492 


dL, Sl, bL 


705, 705, 622 


288, 288, 266 


575, 575, 534 


dR, sr, bR 


676, 676, 667 


276, 276, 276 


550, 550, 549 


fl, T2 


266, 326 


104, 124 


151, 212 


tl, t2 


419, 705 


177, 363 


445, 594 


bi, 62 


620, 670 


264, 278 


533, 550 


Xi ) X2 


239, -468 


53, 257 


192, -329 


=LSP, xlx'i 


126, 239, -457, 464 


27, 65, -220, 263 


102, 191, -316, 324 


Ml, M2, M3=g 


126, 245, 718 


45, 90, 292 


102, 200, 610 


BR{b 57) 


3.52 X 10"* 


2.97 X 10-4 


4.76 X 10-4 


mb{Mx)/mr{Mx) 


0.794 


0.750 


0.794 




0.27 


0.24 


0.22 


B & gaugino purities 


0.99, 0.99 


0.67, 0.87 


0.98, 0.98 



Table 8: Three representative solutions-one with rather light sparticles, and the other two 
with intermediate to heavy sparticles. All masses are given in GeV. Some neutralino and 
chargino masses are quoted as negative. This is merely an indication of the phase of the 
mass eigenstates (expressed as rji by some authors); we include it in case people wish to use 
these number in calculations, but only magnitudes should be considered as experimentally 
relevant, 2 are physical eigenstates, while tL,R correspond to the stop mass-matrix entries. 
Same for the stau and the sbottom. 
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top quark. Similarly, all values of tan (i between one and about 50-60 (perturbative upper 
bound) are still allowed, although the resulting phenomenology often differs considerably in 
the small and large tan [5 regime. On the other hand, significant constraints can be placed 
on the (?TT'i/2) "^o) plane. The region ^ "^o is invariably excluded by requiring either 
electroweak symmetry breaking or neutral LSP, while uiq ^ mi/2 is typically ruled out by 
either EWSB or a lower bound on the age of the Universe. Refs. [|l9t |30|] have argued that 
the region mo ^ rni/2 ~ "i^z (and small tan /? ^ 8) appears to be favored by bounds on the 
proton decay in the simplest SUSY SU{5) model, but we prefer not to rely on an SU (5) 
GUT, so this region is not favored for us. 

We have made a first survey of phenomenological implications for future SUSY searches 
in high-energy experiments. We find reasonable chances for eventually finding a chargino 
at LEP II and the gluino or gauginos at the Tevatron. The light Higgs h has a very good 
chance of being discovered at LEP II but most likely only if its beam energy is pushed close 
to or beyond 200 GeV. On the other hand, the chances are very slim with the currently 
approved ^/s = 178 GeV. The LHC will produce large quantities of all superpartners. 

Several predictions follow from our analysis which could have served to falsify the 
CMSSM before the sweeping supercollider searches for the squarks, Higgs, and the gluino 
are done. 

• We derive a general upper bound on as(m^) < 0.133. For larger mt, and for some 
regions of tan/3 when mt is smaller, there is a lower bound as{mz) ^ 0.117. In the 
regions where the bound is not implied by the physics constraints, it is implied instead 
by the addition of a fine-tuning constraint. GUT-scale threshold corrections may be 
sizeable and modify these limits by several percent. (In this paper we have ignored 
all GUT-scale corrections because we have not yet studied specific unification gauge 
groups and their Higgs structures, although we do assume sin^ 9^{Mx) = 3/8.) 

• Within our parameter space the light Higgs mass is very SM-like and rrih ^ 120 GeV 
for mt = 145 GeV and m/j < 130 GeV for mt = 170 GeV, with somewhat lower values 
usually favored. We also find that mh > 85 GeV for tan /5 > 5. If /i had been 
discovered below about 30 GeV, our entire parameter space would have been excluded. 

• The charged Higgs is always significantly heavier than mw and its discovery should 
not be expected at LEP II and most likely even at the NLC500. Other heavy Higgs 
bosons {H and A) are almost degenerate in mass with H^. If is discovered below 
about 110 GeV, our entire parameter space is excluded. 

• If BR(Z bb) < 0.214 and mt ^ 150 GeV, all solutions are excluded. Similar bounds 
exist for larger mt- 

• The LSP is almost invariably of the gaugino-type (more precisely bino-type), as ad- 
vocated early in Ref. [^ ]. If the cold dark matter is not of this type, almost all 
solutions from this study are exluded. If at least one sfermion (other than the stop 
or sneutrino) had always been lighter than about 80 GeV, there would have been too 
little neutralino DM |110| ]. Furthermore, had the sneutrino been the LSP, the CMSSM 
would not have predicted enough DM. 
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• Several clear patterns and relations among the masses of the Higgs and SUSY particles 
arise which can be tested in futm'e accelerators. 

We have also addressed several related issues recently discussed in the literature. We 
agree with d, I, g |lOl about the need to use two-loop RGEs and to include multiple 
thresholds in considering the running of the gauge couplings. One-loop RGEs and simplified 
one-step approaches can each lead to errors in as{mz) of some 10% (while the experimental 
error is at the 5% level). We emphasize that it is inappropriate to use the so-called effective 
SUSY-breaking scale Mg^sy ™ deriving the GUT mass scale Mx and gauge coupling ax- 
In particular, the value of Mx derived this way can be twice that coming from the full 
two-loop calculation. We have noted that GUT-scale corrections of 10-15% to the relation 
mh{Mx) = rririMx) may have a significant impact on the resulting mass of the bottom 
quark mb{mb) and may allow for this relation to hold for wide ranges of both mt and tan/3. 
We also find it very important to use the one-loop effective Higgs potential particularly in 
the range of large masses ( ^ 1 TeV) where it can lead to qualitatively different conclusions 
than had we simply used the tree level potential. 

Remarkably, we find that for larger mt it is possible to place upper bounds on all 
superpartner masses of 0(1 TeV), without imposing any fine-tuning criterion. For smaller 
mt this is still true for some regions of tan /? near one, but not for all. 

Ultimately the goal is to go from experimental data to a determination of the effective 
Lagrangian of the supersymmetric and (perhaps) unified theory at a scale of 0(10^^ GeV). 
We have shown by example that an effective and perhaps optimal way to do this as data 
becomes available is to systematically reduce the allowed parameter space numerically. Once 
the high-scale Lagrangian is known, perhaps the patterns among its parameters will lead 
toward an understanding of how SUSY is broken and what the underlying theory is. 

Minimal Supersymmetry is a very attractive theoretical framework which makes several 
falsifiable phenomenological and cosmological predictions while at the same time encom- 
passing all the remarkable experimental successes of the Standard Model. The most natural 
ranges for supersymmetric particle masses typically lie above the reach of currently operat- 
ing accelerators (LEP, FNAL) but may be accessible to the upgraded LEP and FNAL, and 
should be successfully explored by the next generation of colliders. 
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Figure 1: Regions in the — tan /? plane consistent with bottom-tau Yukawa unification. 
The region bounded by the sohd hues represents the region of parameter space consistent 
with irih/mT = 1 at Mx for 4.7 < m^°^^ < 5.1 GeV. The region between the dashed hues is 
consistent with mh/nir = 0.9 at Mx- Here we have taken the effective scale of SUSY to be 
90 GeV and as{mz) = 0.120. 

Figure 2: Same as Fig. I but now with as{mz) = 0.112. Notice that the available 
parameter space has increased markedly. 



Figure 3: The running of the sparticle masses from the GUT scale to the electroweak scale, 
for a sample set of input parameters (see "Solution 3" in Table |8[^ater in this paper). The 
bold lines are the three soft gaugino masses, m~, M2 (labelled W) and Mi (labelled B). 

The light solid lines are the squark {qi, qn, ti, ifi) and slepton (Zl, Ir) soft masses, where 
we ignore D-term contributions and the mixing of the stops for this figure. Finally, the 
dashed lines represent the soft Higgs masses, mi and m2 (see Eq. 10), labelled by Hd and 
Hu- The onset of EWSB is signalled by going negative, which is shown on the plot as 
7712 going negative for convenience. 



Figure 4: Plots of the (mi/2, ''t^o) plane showing regions excluded by lack of EWSB (labeled 
E), neutralino not being the LSP (L), the age of the Universe less than 10 billion years (A), 
m+ < 47 GeV (C), BR(6 57) > 5.4 x 10"^ (B), and SM-like lightest Higgs mass 
nifi < 60 GeV (H). We take mt = 145 GeV, sgn^uo = —1, and several representative choices 
of tan/3 and Aq. In window (a) tan/3 = 1.5, Aq/ttiq = 0, in (b) tan/? = 5, Aq/itio = 0, 
in (c) tan/? = 5, Ao/niQ = —2, and in (d) tan/? = 20, Ao/mo = 3. In window (a), the 
regions excluded by each criterion are identified separately, while for windows (b)-(d) only 
the total envelope is shown. For each case, the limit imposed by our fine-tuning constraint 
/ < 50 is shown as a dotted line, disfavoring regions above and to the right of the line. 
Notice the importance of combining several different criteria in constraining the parameter 
space. (Only the most limiting constraints are marked.) Note that in window (a) the 
(mi/2;?7io) allowed region is bounded entirely by the physics constraints, without a fine- 
tuning constraint, though mi/2 extends to larger values than allowed by this constraint (see 
also Fig. 0). 

Figure 5: Same as in Fig. ^ but for mt = 145 GeV, sgn^o = +!> and in (a) tan/? = 5, 
Ao/mo = 0, in (b) tan/? = 5, Ao/mQ = —2, in (c) tan/? = 5, AQ/rriQ = 2, and in (d) 
tan/3 = 10, Ao/mo = -2. 

Figure 6: Same as in Fig. ^ but now with mt = 170 GeV, sgn^uo = —1 and (a) tan/? = 5, 
Ao/mQ = 0, in (b) tan/? = 5, Ao/mo = —2, and in (c) tan/? = 20, A^/mo = 3. In (d) we 
take tan/? = 10, Ao/mo = —2, and sgn^o = +1- Generally, for mt = 170 GeV we find both 
mi/2 and mo bounded from above by physical constraints. 
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Figure 7: The {mi/2,n^o) plane for mt = 145 GeV, tan/3 = 1.5, Ao/mo = 0, and sgn/^o = 
— 1. We show in window (a) the region allowed by all constraints (dark solid) (compare 
Fig. ^). We also mark the bands favored by the CDM (between dashed lines) and MDM 
scenarios (between dotted lines) and the limit imposed by our fine-tuning constraint / < 50 
(light solid). We also plot in window (b) as{mz) at 0.120 through 0.132 in steps of 0.002, 
decreasing for larger {mi/2,iT^o) with 0.120 solid, 0.126 dashed and any others dotted; in (c) 
m ± (solid), m~ (dashes), and m~ (dots) at 45, 80, 150, and 250 GeV. (In the two latter 
cases only the last three values occur in the graph.) In window (d) we plot uiy. (solid) at 18 
(thick), 45, 75, 100, 125, and 150 GeV. We also display the gaugino purity {Zf^ + Zf^, dots) 
of 0.8 and 0.9 increasing to the right. In window (e) we plot m~ (dots) and average niq 
(other than m^) (solid) at 250, 500, 750, and 1000 GeV, and in window (f) m/j between 60 
and 130 GeV in 5 GeV intervals (60, 120 dark solid, 90 dark dashes, all others light solid). 
Mass contours in each window increase with increasing (mx/2i"^o)- 

Figure 8: Same as in Fig. |^ but for nit = 145 GeV, tan/3 = 5, Ao/mo = —1, and 
sgn//o = — 1. Window (a) also demonstrates the importance of including final states other 



than // (in this case hh) in the calculation of ^^^^o (^^^ Section 8.1.7). 



Figure 9: Same as in Fig. ^ but for mt = 170 GeV, tan/3 = 20, Aq/ttiq = 1, and sgn/^o = 
-1. 

Figure 10: Region in the (mi/2 5'"T'o) plane ruled out by the LEP bound on the sneutrino 
mass, mi) > 43 GeV, for both small and large tan/3 and by the CDF bound m~ > 141 GeV 
(dashes, cascade decays neglected). The chargino mass bound m^± > 47 GeV often (but 
not always) leads to additional excluded regions. 

Figure 11: For the case presented in Fig. |8| {mt = 145 GeV, tan/? = 5, Ao/mo = —1, and 
sgn//o = —1) we delineate the regions close to the Z (window (a)) and h (window (b)) poles 
in the process XX ~^ ff where our calculation of O^^/iq cannot be trusted. (The effect of 
other poles is much less significant.) In window (a) \m-^ — mz/2\ = 25 (dots), 15 (dashes) 
and 5 GeV (solid). For \m^ — mz/2\ ^ 15 GeV our calculation of O^/ig is sufficiently reliable. 
In window (b) the same for the light Higgs h. 

Figure 12: Scaling behavior of the fine-tuning constant with the SUSY scale. The solid line 
represents a fine-tuning of 50 which typically corresponds to m~, m~ ^ 1 TeV. The other 
lines are (left to right) for 1, 10, 100, 500, and 1000. Here we have taken mt = 145 GeV, 
tan/3 = 5, A^/mo = — 1, and fi < 0. 

Figure 13: Scatter plot of (a) m~, (b) m-^, (c) mo, and (d) mi/2 vs. fine-tuning for 
solutions consistent with all applied constraints. Notice that the cut / < 50 typically gives 
sparticle masses m~, m~ < 1 TeV but in some cases (all of which have large tan /3) they can 
be significantly heavier. 
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Figure 14: Scatter plot of Bo/niQ vs. Ao/niQ for all allowed solutions (COMPASS) with 
nit = 145 GeV. The quantized appearance is due to numerical sampling and is not signifi- 
cant. 



Figure 15: Plots of the (W'i/2;^o) parameter space showing regions excluded by lack of 
EWSB (labelled E), LSP not being the neutralino (L), and the age of the Universe (A), 
for (a) no 1-loop contributions to Vyiggs and (b) leading 1-loop contributions to ^Higgs- See 
text for discussion of full 1-loop contributions to VHiggs. For these plots we have taken 
rrit = 145 GeV, tan/? = 5, Ao/niQ = — 1 and sgn/^o = — 1 (compare Fig. Ila). 



Figure 16: Lowest rrih versus .^rrTilnf^ for all acceptable solutions with mt = 145 GeV. 
We allow rrih < 60 GeV for the purposes of this graph only. 



Figure 17: Plot of nih vs. tan /? for all acceptable solutions with mt = 145 GeV. The 
solid vertical bands express the range in rrih for a given tan (3. The dotted lines show the 
clear envelope of rufi vs. tan /3 that we obtain in the CMSSM. Note that if mh ^ 85 GeV 
then tan /3 ^ 5. The discretization of tan (3 is merely from numerical sampling and is not 
physically significant. 



Figure 18: Scatter plot of nih vs. \/ \Bfj,\ for all acceptable solutions with nit = 145 GeV 
and tan/3=5. Note that ^/JBJ^\ is usually larger than mz- 
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Figure 19: Histogram of rrih for all acceptable solutions with mt = 145 GeV. 



Figure 20: Percentage of acceptable solutions for rrit = 145 GeV with h being detectable 
at LEP II versus center of mass energy (GeV). 



Figure 21: Scatter plot of tan/3 vs. niA for all acceptable solutions with mf=145 GeV. The 
band structure is due to numerical sampling and is not physically significant. The dotted 
triangular region is the approximate region in which it is difficult to detect at least one 
Higgs jlOOl ]. 



59 



Figure 22: Scatter plot of nih vs. for all acceptable solutions with mj=145 GeV. The 
various bands correspond to different choices of tan /?. The gap in the lower portion of the 
plot is due to solutions with 1.5 < tan/? < 3 which were missed due to our finite grid. 



Figure 23: The ratio of the top Yukawa to the bottom Yukawa couplings at the GUT scale 
vs. tan/? for all acceptable solutions with mt = 145 GeV. We see that (^7^) tan/? ~ 60 is 
not a bad description, so the running cannot help much to account for the large size of the 
mt/rrih ratio; it must be imposed either via or via tan/?. If tan/? is not large, as we 

might expect, then a large (fj^) is required as input. 

Figure 24: Histogram of BR(5 — > 57) for all otherwise acceptable solutions with mt = 
145 GeV. The dotted lines indicate the upper and lower bounds on BR(6 57) imposed 
as discussed in the text. The SM prediction is 3.1 x 10~^ corresponding to the formulas of 



Ref. |7|. More recent QCD correction estimates @ lead to a SM BR(6 sj) ~ 4.3 x 10"^ 



Figure 25: Scatter plot of BR(6 — > 57) vs. mjj± for all acceptable solutions with mt = 
145 GeV. The faint banding visible in the figure is from numerical sampling and is not 
physically significant. 

Figure 26: BR(6 57) vs. tan/? for all acceptable solutions with mt = 145 GeV. In 
order to demonstrate the density of points, tan /? is slightly smeared around its numerically 
sampled value. 

Figure 27: Regions of the {mQ,mi/2) plane to be explored by SUSY searches at LEP II 
if mt = 145 GeV. Crosses (empty boxes) mark the regions that will always (never) be 
accessible to LEP II for any combination of input parameters. Filled boxes mark the regions 
accessible for some combinations of parameters. Empty regions are excluded by our set of 
constraints. In window (a) we show the combination of possible (direct) SUSY searches at 
LEP II as described in the text. In window (b) we show the LEP II search potential for 
the chargino Xi^ alone, and in window (c) and (d) for the lightest Higgs assuming ability 
to find h with masses m^ < 80 GeV and 110 GeV, respectively, {h is a very SM-like Higgs.) 
Very large values of mi/2 are disfavored by the fine-tuning constraint. 

Figure 28: Tn + — m^o vs. m + for all acceptable solutions with mt = 145 GeV and 

Xl Al Xi 

m+ < 120 GeV. The line represents m^± = 2m^o which is approximately true for a 

Xl Xl Xl 

gaugino-like LSP. While the line is an approximate description of the results it is not accurate 
enough for detailed use. 

Figure 29: Histogram of lepton or jet energy from decay for acceptable solutions with 
m + < 120 GeV. 

Xl 
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Figure 30: Histogram of Rb = r^g/Fhad for all acceptable solutions with rrit = 145 GeV and 
0.16 < ^x^o < 0.33. The central dotted line is the measured value of Rb and the outside 
dotted lines are the one standard deviation errors on the measurement. The dashed line is 
the Standard Model calculated value of R^ given rrit = 145 GeV. 



Figure 31: Scatter plot in the plane (/x, M2) for the acceptable solutions with mt = 145 GeV. 
Notice a large concentration of points below the diagonals M2 = which shows that the 
LSP is gaugino-like in most of the solutions. 



Figure 32: Scatter plot of the LSP wavefunction for the acceptable solutions for which 
0.16 < Q^^l < 0.33 (MDM scenario), and for m* = 145 GeV and tan/? = 10. Each solution 
contributes scatter plot points corresponding to its B, W, Hd, and components. Note 
that the LSP is mainly B in all these solutions. 



Figure 33: Histogram of ^x^o acceptable solutions (with an initial selecting cut 

^x^o < 1 imposed) for mt = 145 GeV. Notice a strong peak around « 0.1 showing 

that in the CMSSM the MDM scenario is somewhat more favored relative to the CDM one. 



Figure 34: ^x^o "^^^ (^) "^O) (b) 1^1/2, (c) tan/3, and (d) fi, for otherwise acceptable 
solutions with mt = 145 GeV. The bound f^x^o < ^ comes from the age of the Universe of 
10 billion years or more. The ranges 0.16 < ^^x^o ^ ^-^^ ^^'^ ^-^^ ~ ^x^o ~ ^-^ favored 
by the MDM and CDM scenarios, respectively. The banding in mo, mi/2 and tan/3 is from 
numerical sampling and is not physically significant. 



Figure 35: Scatter plot of mass vs. particle type for all acceptable solutions in our data 
set with mt = 145 GeV and / < 50. The horizontal bands are due to numerical sampling 
and are not of significance. 
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